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ABSTRACT  N^rious  models  of  Navy  walc;froiit  structures  arc  presented  and 
employed  to  assess  their  efTectiveness  in  predicting  the  vulncrahilily  of  these  systems  to 
opcratiorul  and  earthquake  loads.  The  emphasis  is  on  application  of  three-dimensional 
finite  element  models  to  reinforced  coiKreic  structural  analysis.  Results  indicate  that 
dynamic  arulysis  is  important  in  the  calculation  of  the  structura.'  response  of  both  a 
drydock  and  a  blocked  vessel  in  drydock.  The  primary  nahiral  frequencies  and  natural 
mode  shapes  for  a  drydock  arc  longitudinal  deformation  modes  which  cannot  be  pre¬ 
dicted  by  current  methodology  based  upon  statically  equivalent  analysis  of  two-dimen¬ 
sional  models.  Further  results  from  nonlinear  analysis  provide  new  insight  into  the 
behavior  of  the  drydock/caisson  seal  for  hydrostatic  loads,  w  hich  is  dramatically  differ¬ 
ent  fnjm  that  which  underlies  current  design  and  maintenance  procedures.  Similarly, 
results  from  a  three-dimensional  nonlinear  static  analysis  of  a  .scale  model  of  a  rein¬ 
forced  coiK'ietc  pier  deck,  subjected  to  punching  shear  failure  loads,  arc  shown  to 
comp.*e  well  with  cspcimcntal  data.  MfxJer.i  three-dimensional  finite  eleriKiit  tech¬ 
nology  is  appropriate  foi  analysis  of  waterfront  stmcturcs. 
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1*.  Dl9TRIRUTK>MCOOC 

IX  AMTHACT  (MA>>nun  aoo  •»<*; 

Various  models  of  Navy  waterfront  structures  arc  presented  and  employed  to  assess  their  effectiveness  in  predict¬ 
ing  the  vc'ncrability  of  these  systems  to  operaticnal  and  earthquake  loads.  The  emohasis  is  on  applicafio;.  of  three- 
dimensional  nhitc  elcriKnt  models  to  reinforced  concrete  structural  analysis.  Results  indicate  that  dynamic  analysis 
is  important  in  the  calculation  of  the  structural  response  of  both  a  drydock  and  a  blocked  vessel  in  drydock.  The 
primary  natural  frequencies  and  natural  mode  shapes  for  a  drydock  are  longitudinal  deformation  modes  which  cannot 
he  predicted  by  current  methodology  based  upon  statically  equivalent  analysis  of  two-dimensional  model.s.  Further 
results  from  nonlinear  analysis  provide  new  insight  into  the  behavior  of  the  dry  dock/caisson  seal  for  hydrostatic 
loads,  which  is  dramatically  differcni  from  that  which  underlies  current  design  and  maintenance  procedures.  Simi¬ 
larly,  results  from  a  three-dimensional  nonlinear  static  analysis  of  a  scale  model  of  a  reinforced  concrete  pi'^r  deck, 
subjected  to  punching  shear  failure  loads,  are  shown  to  compare  well  with  experimental  data.  Modern  three- 
dimensional  tinitc  element  technology  is  appropriate  for  analysis  of  waterfront  structures. 

14.  MMUtCTTEmS 
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reinforced  concrete  pier  deck,  concrete  material  models,  waterfront  facilities 
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aiibnt  Ant  la  1  capital  tnvAatmrnt.  and  a  anhatatitlal  atratcftic  raaource  for 
tha  Navy.  Thay  nay  alao  ba  claanlflad  aa  aRaantial  and/or  hlgh-rlak 
facllltlaa.  Corraapondlngly ,  thay  rnqiilra  apar.lal  attantlon  In 
anginaaring  daalgn,  har.ard  mitigation  and  matntananca  maaatiraa,  and 
vnlnarabil Ity  atudlaa.  At  tha  aama  lima,  anbstantlal  tachnology  In 
structural  anatyaia  capability  exiata  dua  to  advances  In  computational 
structural  machanlcs.  Unfortunately,  this  existing  tachnology  la 
ralativaly  untapped  regarding  approaches  to  design,  analysis,  and 
vulnerability  asaassment  of  the  Navy's  facilities.  Experience  In  Its 
application  must  be  gained.  Tills  problem  l.s  adilresaed  in  this  report. 


OBJECTIVE 


The  objective  of  this  study  was  to  demonstrate  the  effectiveness  of 
existing  modern  computer-based  structural  analysis  methods  by  applying 
them  to  several  current  technical  problems  involving  naval  waterfront 
facilities. 


BACKGROUND 

Uninterrupted  operation  of  drydocks  and  piers  Is  essential  to 
readiness  of  fleet  combatants.  During  war,  waterfront  systems  support¬ 
ing  repair  of  fleet  combatants  would  be  threatened  by  air  forces  and 
underwater  swimmers.  Global  Wargame  B6  Indicated  that  fleet  operational 
limitations  could  result  from  current,  graving  drydock  vulnerability. 
During  peacetime,  ship  repairs  can  be  catastrophically  Interrupted  by 
damage  incurred  from  a  major  seismic  event.  The.re  Is  a  5  percent  annual 
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probfihility  of  a  mfljor  nftlsmir.  o.vont  {>7  R(r.ht**r  scfllo)  with  tho  potential 
of  severely  demaglng  graving  drydor.ks  at  Puget  Soiinfl,  Mare  Island,  Pearl 
Harbor,  Hunters  Point,  and  I,ong  Bearh  Naval  Shipyards,  and  San  Diego 
Naval  Station.  The  minimum  replacement  <:f>sl  for  a  small  drydock  has 
been  estimated  to  be  $140  million.  The  potential  damage  to  docked  ships 
and  the  impact  of  disruption  of  fleet  overhaul  cycles  Is  Incalculable. 

For  the  sake  of  brevity,  this  report  excludes  .consideration  for  the 
Important  problem  of  Naval  SI. ore  Facilities  subject  to  blast  loads. 

This  Is  the  subject  of  a  separate  report  entitled,  "Three-Dimensional 
Structural  Analysis  Methodology  for  Navy  Explosive  Safety  Facilities  -  A 
Technology  Assessment"  (Shiigar,  et  al.,  1002). 

STRUCTURAr.  ANALYSIS  SUBSYSTEM  HODEI-S  STUDIED 

One-dimensional  continuous  models  of  drydocks  and  their  natural 
frequencies  and  mode  shapes  are  discussed  and  demonstrated  initially  in 
this  report.  Studies  of  these  preliminary  models  provide  orientation 
and  perspective  on  the  behavior  attainable  from  more  complex  three- 
dimensional  finite  element  models. 

An  assessment  is  made  of  two-dimensional  models  of  drydocks.  These 
models  may  be  said. to  be  the  present  technical  basis  for  structural 
engineering  aspects  of  drydock  certification.  Tliey  are  useful  in  pro¬ 
viding  needed  experience  in  application  of  finite  element  technology. 
However,  the  assumptions  upon  which  they  are  based  are  ofte?i  not  satis¬ 
fied. 

The  safety  of  a  blocked  vessel  in  drydock  isthe  paramount  concern 
of  the,  drydock  certification  process.  It  follows  that,  technology  per-, 
mitting,  the  blocked  vessel  be  viewed  as  a  flrydock  subsystem  so  that  its 
dynamic  response  Is  coupled  with  the  dynamic  response  of  the  drydock 
itself  allowing  the  response  to  be  calcu-Iated  more  accurately;  A 
simplified  two-dimensional  lumped  parameter  model  of  a  blocked  vessel  in 
drydock  is  presented  and  its  dynamic  response  is  computed  relative  to  a 
simplified  strong  motion  earthipiake  load  to  demonstrate  this  modeling 
concept .  . 
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To  dftMonitmtd  th*ilr  rtnnlytfci*!  pf  r<»rt  thrpp'dlwprm lon*l 

flnlt«  •Iniupnt  mod«l«  of  drydocks  /ind  drydo*;k  fiMb«y«t pmsi ,  particularly  a 
■  tppl  calaaon,  are  inveat  (gated .  The  effective  cotiat  met  li>n  of  three- 
dlieeneiona  finite  element  mcadela  la  dlacuaaed. 

A  del  ailed  three-dlmetia Iona  1  finite  elemenl  model  of  a  drydock 
calaaon  auhayntem  waa  developed.  Analyses  using  this  model  are  pre« 
Rented.  Tliey  Include  a  study  of  the  contact  forces  to  which  the  seal  la 
subjected  due  to  typical  hydrostatic  load  condlt  ion.s.  Alao,  natural 
vibrations  of  the  calaaon  model  are  presented  and  dlsciiased. 

A  simplified  three-dimensional  finite  element  model  of  a  reinforced  , 
concrete  drydock  subsystem  was  constructed  based  uf«)n  plate  bending  finite 
elements.  The  model  accounts  for  the  varlaltle  thickness  of  the  drydock's 
sidewalls  and  endwalls  and  alao  Includes  the  steel  caisson  as  part  of 
the  subsystem.  In  a  study  of  the  drydock  subsystem's  natural  vibration 
characteristics,  the  first  ten  natural  f refjnenc ies  and  natural  mode  shapes 
were  calculated  to  demonstrate  that  the  basic  frequencies  of  response 
and  modes  of  deformation  likely  to  participate  in  any  dynamic  response 
to  earthquake  excitation  are  different  from  those  predicted  by  two- 
dimensional  models. 

Other  waterfront  facilities  that  impact  fleet  readiness  are  piers 
and  wharfs,  A  three-dimensional  nonlinear  analysis  of  a  1/3-scale 
laboratory  model  of  a  reinforced  concrete  pier  was  conducted.’  Tills 
analysis  evaluates  the  ability  of  current  finite  element  technology  in 
accurately  reproducing  the  behavior  of  highly  nonlinear  materials  such 
as  reinforced  concrete.  To  this  end,  concrete  .and  steel  rebars  were 
modeled  separately  in  what  is  known  as  a  discrete  reinforcement  model. 

The  development  of  cracking  in  the  concrete  was  followed,  together  with 
the  associated  stiffness  degradation  and  nonl inear  load-displacement 
history. 


STRUCTURAL  ANALYSIS  SOFTWARE 

Though  the  demonstration  analyses  de.scribed  above  are  the  primary 
subject  of  this  report,  the.y  are  reg.irded  as  preliminary  to  the  development 
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of  a  specialized  nonlinear  atructisral  analyals  software  nystem  for 
selected  large  striictiiral/geotechnlcal  syatema  in  naval  facilities. 

The  overall  goal  la  to  adapt  appropriate  modern  methods  In 
computational  mechanics  and  combine  them  with  several  more  recent 
products  from  Navy-sponsored  basic  research  In  structural  modeling  to 
produce  a  special  purpose,  advanced  nonlinear  structural  analysis 
system.  Some  of  these  basic  research  products  and  their  transition 
potential  were  reviewed  In  "An  Evaluation  of  Numerical  Algorithms  for 
the  Nonlinear  Dynamic  Analysis  of  Large  Soil  Structure  Systems"  (Bayo, 
1987). 

The  flexibility  to  effectively  and  'exped 1 t Ions ly  address  and 
determine  the  essential  structural  behavior  of  any  critical  facility 
subjected  to  strong  motion  earthquake  and  severe  blast  loads  is  a 
desirable  goal,  and  the  proposed  special  purpose  software  system  will 
facilitate  achievement  of  that:  goal. 

While  commorclally-avallable  technology  is  potent.  It  very  often 
does  not  strictly  apply  to  complex  Navy  systems.  As  a  result,  the 
technical  problem  must  Inevitably  be  modified  to  accommodate  the  problem 
specification  requirements  of  existing  commercial  structural  analysis 
softw'are  products.  These  products  have  been  developed  for  a  very  large 
market  place  In  which  the  Navy  is  but  a  single  customer.  They  provide  a 
general  purpose  capability  and  hence,  they  often  do  not  strictly  apply 
to  specific  Navy  problems.  Experience  has  born  this  out.  The  contracting 
cycle  is  not  suited  to  expeditious  development  and  procurement  of  a  new 
or  modified  commercial  capability.  Hence,  timely  acquisition  of  the 
required  capability  most  often  does  not  occur.  It  is  also  true  that  lii 
nonlinear  problems  the  required  capability  often  cannot  be  well-defined 
until  several  attempts  to  solve  the  f)’"obl,em  have  been  made.  All  this 
pKDlnts  to  the  need  for  enhanced  response  via  a  Flexible  capability  in 
nonlinear  structural  analysis.  A  software  framework  designed  to  provide 
gi’loance  toward  achieving  this  c.cpability  was  developed.  Some  of  the 
Issues  considered  are  briefly  summarized  in  the  following  paragraphs.  A 
full  description  is  presented  in  "A  ooftware  Development  Specification 
for  Nonlinear  Structural  Analysis"  (Landers,  1990). 
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Structural  anglnaarn  ah«til*l  •xptolt  ri»«  •nt  a<lv«nc«»  In  coaiputar 
acianca  with  ragard  to  hardware  and  software.  I^IxiatinK  batch  ortantad 
•nvJ romaanta  ar«  conduclva  to  cr»at)fi|{  <r.advart.»nt  arrora  In  Input  data 
which  nagala  long  coatly  analyria  coaipntar  rnna .  Furthar,  ihay  maka  tha 
davatoiiaiant  and  taatlng  oT  iiaw  mathodologlaa  and  algorlthma  vary  diffi¬ 
cult,  Kvan  wl  th„  1  nt aract  1  va  taxt  adltora,  tha  <  f>nat  rtjct  fon  and  debugging 
of  naw  concapta  and  idaaa  oftan  fnvf>lvaa  atany  itaratlona  of  tha  ”adlt- 
coatplla-dabug”  eye  la, 

Low-coat  work  atatlona  aquippad  with  powerful  12-blt  procaaaora, 
hlgh-raaolutlon  graphic  diaplaya,  and  Inaxpanalva  networking  facilltiaa 
are  appearing  rapidly  while  the  coat  of  raw  procosaing  power  la  dropping 
rapidly.  Further,  and  a ignl f leant ly ,  vendora  are  h»  ^Innltg  to  coma  to  a 
conaansua  on  atandarda  for  the  toola  they  provide  for  Intermachlna 
communication,  gr.iphlca,  it'd  data  baaes. 

Unfortunately,  advanc^c;  in  atri'cturil  engineering  aoftwara  develop¬ 
ment  have  not  kept  pace  with  the  rapid  changes  in  tiie  computer  market 
p’ace.  Some  attempts  are  documented  wherein  existing,  we  1 1 -respected, 
awnollthlc  finite  element  software  systems, traditional ly  run  on  large 
computers  have  been  "down  loaded"  to  workstation  environments.  However, 
little  Innovative  software  is  available  to  the  research  and  development 
community.  This  group  requires  a  compiitat ional  environment  that  is 
re.sponsive  to  changes  in  the  state  of  the  art  in  computational  mechanics 
and  software  development.  Since  new  rr.'-.niques  and  algorithms  must  be 
tested  and  debugged  as  they  are  implem-nted,  a  systemthat  provides  a 
high  degree  of  flexibility  and  interaction  is  required. 

A  description  of  a  software  environment  for  application  to  struc¬ 
tural  engineering  is  provided  In  the  i.ander.s  (I'l'lO)  report.  The  system 
is  primarily  Intended  for  the  design  and  implementation  of  new  methodolo¬ 
gies  and  techniques  in  finite  element  analysl.s,  but  it  can  be  extended 
to  production  situations  as  well. 
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DCHONSTRATION  OF  STRUCTURAI.  ANALYSIS  NODEI^S 


ONE-DINRNSIONAL  CONTINUOUS  BF.AN  HOOF.I.S  OF  DRYDOCRS 

Since  A  drydock  h’AR  a  long  Alender  Ahepe,  (t  Ja  tempting  to  treAt 
It  AS  A  continiioiiA  beAm  for  purpOAeA  of  a  prellminAry  AtructurAl  dynAmics 
study.  Indeed  eArly  design  And  AnstyAis  spproAclieA  for  drydocks  were 
bssed  on  elementnry  besm  lormules  (URS/John  A.  Bliime  and  Associates, 

1978;  Woodward-Clyde  Consultant.s ,  1979;  and  Moffat  ard  Nlchol,  Engineers, 
1981).  These  early  approaches  were,  however,  static  analyses  where  the 
drydock  walls  were  treated  as  vertical  cantilever  l-eams.  It  Is  surprising 
that  certain  fundamental,  readily  available  results  for  one*dlmensional 
contlniipuA  beam  natural  frequencies  and  natural  vlbra<  Ion  mode  .shap.^s 
were  not  considered  In  these  early  studies.  In  this  section,  these  data 
are  presented  to  demonstrate  their  effectiveness  in  two  respects.  First, 
the  longitudinal  mode  shapes  of  drydccks  can  be  estimated  by  one-dimensional 
continuous  beam  models.  Second,  these  data  serve  to  guide  the  development 
of  more  advanced  analyses  based  on  three-dimensional  finite  element  models 
of  drydocks. 

Three  separate  continuous  beam  models  are  presented  and  discussed, 
each  Increasing  In  complexity  and  scope.  The  data  presented  fer  the 
first  and  second  models  are  essentially  available  In  various  reference 
books  on  structural  dynamics  such  as  those  listed  in  the  bibliography 
provided  In  Appendix  A.  The  third  model  is  a  newly  derived  modification 
of  the  other  models  to  account  for  the  translational  Inertia  of  the  drydock 
endwall  and  drydock  caisson. 

Properties  for  the  one-dimensional  continuous  beam  models  are  cal¬ 
culated  based  upon  an  Idealized  cross  section  of  drydock  No.  6  at  Puget 
Sound  Naval  Shipyard  (NSY).  This  particular  drydock  was  selected  because 
it  Is  a  relatively  new  drydock  with  a  thin  wall  and  a  thin  floor  structural 
design.  Its  engineering  design  and  construction  is  well  documented  in 
the  open  literature  by  Zola  and  Boothe  (1960).  The  idealized  cross  section 
and  associated  data  used  are  presented  in  Figure  1.  The  material  model 
is  linear  elastic  and  the  mass  density,  and  modulus  of  elasticity, 
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h«650in. 

w  ■  2,250  In.  PcA  ■  84.8  lb  •  t^/ln.^ 

t2  -  70  In. 

tj  .  120  In.  Egly  «  5.68  X  10^8  |b.|n.2 

t^«110ln. 


Figure  1 .  idealized  Cross-Section  Properties  of  Orydock  6  at  Puget  Sound  NSY 


E^,  are  nominal  values  for  concrete.  The  product  of  and  the  cross- 
sectional  area  A  represents  the  mass  per  unit  length  of  the  drydock,  and 
the  product  of  and  the  second  montenl  of  area  about  the  neutral  axis 
represents  the  longitudinal  bending  rigidity  of  the  drydock. 

Free-Ended  Beam  Model 

In  the  free-ended  beam  model,  the  drydock  Is  assumed  to  have 
free-end  conditions.  Thi.s  is  also  sometimes  referred  to  as  a  "flying 
beam"  model.  The  effects  on  natural  frequencies  and  mode  shapes  of  the 
drydock  due  to  the  surrounding  soil,  and  foiind.al  ion  are  Ignored  In  this 
model . 

A  tabulation  of  natural  frequencies  and  mode  shapes  for  this  model 
and  other  continuous  beams  can  be  found  in  the  comprehensive  reference 
by  Blevins  (1979).  The  results  for  the  present  drydock  model  are 

presented  in  Figure  2.  Here  I.  Is  the  length  of  the  drydock  and  f  is 

th  ^ 

the  natural  frequency  for  the  n  natural  mode  ot  ViUiatlon.  The 

formula  for  this  frequency  is  given  in  the  figure.  The  first  five 
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non*«»ro  nutuntl  rtf*»  fili»o  ilpp<rt»»fi.  Th#««5«  r«>!tuttN  naturally 

pertain  to  tonaf h^nflfi’g  In  tlio  <}ry«<n«:k  siHpwalla  and  floor  (aa 
cnntra-ftffd  with  vertical  bonding  in  t  ho  walla  from  a  cant  i  lavorad  wal  1 
modal)  Tha  dynamic  rraponaa  for  d isp I ar amant a  and  atraanaa  In  tha  dry* 
dock  ralativa  to  a  pra«  :rlh«‘d  dynamir  load  can  ha  calculatad  lining  a 
I  inaar  combination  of  thana  fnndamantal  daformation  nhapaa  accor''ing  to 
tha  mrvlal  aiiparpoa I t ion  mathod  in  atrnrfnral  dynamica  thaory  an  daacrlbad 
by,  fpr  aaampla,  flloiigli  and  Panzian  (l*)?^).  Tha  particular  combination 
of  mrvta  Mod  tha  magnltiida  of  tha  taajronRa  dapand  on  tha  fraquancy 

contant  and  spatial  dist r ibnt ion  of  (ha  piaRcribad  dynamic  load. 

In  addition  to  naglartlng  tha  affact  of  tha  .•soil  and  foundation, 
thaaa  raaiilts  ara  basad  on  tha  assumption  that  vibration  of  tha  drydock 
is  syrnmat  r  i  ca  1  about  a  vartic.al  plana  ihtough  tha  drydock  longitudinal 
r.antarlina,  Consaquant  1  y ,  othar  vibr.ation  modas  such  as  torsional 
daformation  Hwidas  (longitudinal  twisfing  about  Ilia  c.antar  1  ina)  and  bend¬ 
ing  sKvias  in  tha  horizontal  plana  (about  tba  strong  axis)  are  pra.cludad 
from  thasa  data.  To  ratain  thas'a  mo«las,  thraa-d  imans  lonnl  models  must 
ba  considarad.  Howavar.tha  affact  of  an  assiimad  elast Ic  foundation  for 
tha  drydock  is  aasily  investigatad  by  ona-diman.s  ional  continuous  beam 
sndels  as  shown  In  the  next  modal. 

El  ant  ica  My -Supported  Free-F.ndad  gaam  Model 

In  the  alastical ly-siipportad  fraa-andad  baam  model,  the  drydock  is 
assiimad  to  ha  continuously  supportad  along  its  langth  by  an  elastic 
foundai  ion  as  shovTi  in  Figure  3.  Tlia  modulus  of  tha  foundation,  E^.,  has 
units  of  force  per  unit  area,  and  tha  .stiffness  mo<hiIiis,  k^,  has  units 
of  force  par  unit  langth.  Tha  result ing  natural  frequencies  and  natural 
mfwias  for  this  system  are  shown.  Tliasa  .lata  are  also  taken  from  tabula¬ 
tions  in  Blevins  ( 1979) . 

The  natural  fraquancfas  of  the  system  increase  with  the  square  root 
of  the  foundation  modulus,  F.^,  as  shown  by  the  formula  given  in  Figure  3 
(note  that  this  is,  the  same  as  the  formula  in  Figure  2  when  E^  is  zero). 
However,  the  natural  mode  shapes  of  the  continuous  beam- foundation  model 
are  identical  with  the  natural  mode  shapes  of  the  continuous  beam  model 
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Frequencies  and  Modes  for  a  Free-Ended  Beam  on  an  Hastic  Foundation. 


in  the  absence  of  the  foundation.  ThI.s  well-known  result  was  first 
published  by  Stafford  (  1967'*.  Simply  stated,  the  shapes  of  the  natural 
longitudinal  bending  mod js  of  the  drydock  are  Independent  of  the  founda¬ 
tion  modulus,  ar cording  to  this  model. 

ElasLlcal  ly-Supported  Fre.c-Ended  Ream  With  Tip  Masses 

To  account  for  the  trans lat Iona  1  Inertia  of  the  drydock  endwall  and 
the  drydock  caisson,  the  previous  continuous  beam  model  was  modified  to 
include  concentrated  masses  at  either  end  of  the  continuous  beam,  m^ 
and  m^,  respectively.  This  model  is  shown  in  Figure  4. 

From  the  previous  model,  it  was  seen  that  the  presence  of  the 
elastic  foundation  caused  an  increase  in  the  natural  frequencies  of  the 
system  in  a  very  simple  way.  This  is  also  true  of  the  present  model. 

It  can  be  shown  (see  Stafford,  1967)  that  the  natural  circular  frequency 
$2,  of  a  beam  on  an  elastic  foundation  can  be  determined  by  taking  the 
square  root  of  the  sum  of  the  squares  of  the  "rigid  body"  natural  cir¬ 
cular  frequency,  given  by^JE^/(p^A) ,  «..d  the  natural  circular  frequency 
of  the  beam  in  the  absence  of  the  elastic  foundation,  «.  The  relation 
is  general  and  Independent  of  boundary  conditions  and  holds  for  any 
natural  mode  of  vibration  n.  Thus, 


There  are  no  published  tabulated  natural  frequency  data  for  the 
continuous  beam  model  with  tip  mas.ses.  Therefore,  the  characteristic, 
equation  governing  the  natural  frequencies  and  mode  shapes  and  its 
solution  are  presented.  The  theory  used  to  derive  this  equation  is 
explicated  well  in  Craig  (1981),  for  example.  The  resulting  charac¬ 
teristic  equation  is; 

2  tanh  XL  tan  XL 

+  (Mj  +  M^)  X  (tan  XL  -  tanh  XL) 

+  - - - =  1 

cosh  XL  cos  XL 
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y.  v(x,  t) 


Rgure  4.  Elastically  Supported  Free-Ended  Beam  Model 
with  Tip  Masses:  1-D  Drydock  Model. 


where:  = 


and  = 


^2 


and 


x"  = 


p  Aw 

E  I 
c  y 


The  tip  muss  ratios,  and  M2,  are  tip  masses  normalized  on  the 
mass  per  unit  length  of  the  beam.  It  is  noted  that  when  the  tip  mass 
ratios  are  zero,  the  characteristic  equation  for  the  "flying  beam"  is 
recovered  from  this  equation.  This  equation  is  (see  Rogers  1959,  for 
example): 


cosh  XL  cos  XL  =  1 

The  "flying  beam"  natural  frequencies  (indicated  in  Figure  2)  are 
determined  from  the  roots  of  this  equation. 

It  is  not  the  tip  masses  m^^  and  m2,  but  the  tip  mass  ratios  and 
M2,  that  matter  when  determining  the  natural  circular  frequencies  for 
the  longitudinal  bending  modes  of  the  drydock. 

For  purposes  of  solving  the  characteristic  equation  for  its  roots, 
let  a  =  XL.  Then,  the  characteristic  equation  becomes: 
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a  tanh  (a)  tan  (a)  + 


a  (tan  a  -  tanh  a) 


cosh  (a)  cos  (a) 

For  the  case  of  drydock  No..  6  at  Pngot  Sound  NSY,  the  above  dimen.s  ion  less 
coefficients  are  evaluated  as: 


^1  ^2\  -12 
I - i =  5.1176  X  10 

— 2|  ^  4.1270  X  10  ® 

The  solution,  by  numerical  iteration,  of  the  characteristic  equation 
for  the  first  root  gives: 


=  4.57565 


or 


'1  = 


4.57565 


Therefore,  the  first  (nonzero)  circular  natural  frequency  is: 


"1 


E  I 


=  6.5524  rad/s 

)  A  1 


The  graph  of  the  iterative  solution  flepicting  the  root  a^^  is  shown  in 
Figure  5.  The  corresponding  value  when  the  tip  masses  are  zero  is  (from 
Figure  2):  ' 


6. 916  rad/s 


which  is  higher  than  the  value  when  tip  masses  are  included.  Thus,  inclu¬ 
sion  of  the  tip  masses  in  the  continuous  beam  model  causes  a  6  percent 
decrease  in  the  fundamental  natural  frequency. 
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y^{a)  ■  5.1176  X  1(r'*2  tanh(a)  tan(a) 

+  4.1270x10*^a(tana-tanha)  ♦  - ^ - 

cos)i(a)  co8(a) 

VjW  « 1 


Using  the  lower  VHlue  of  the  one-dimensional  continuous  beam 
model  estimate  of  the  first  natural  fr«M|ueiicy  of  drydock  No.  6  at  Puget 
Sound  NSY  is: 

fj  =  Wj/(2Tr)  =  1.04  Hz 

This  frequency  is  indeed  clone  to  the  range  of  frequencies  containing 
substantial  energy  in  typical  earthquake,  ground  motion  records.  However, 
this  value  Ignores  the  effect  of  the  soil  foundation. 

Various  values  of  the  soil  modulus,  E^,,  were  compiled  from  Wu  (1967) 
to  demonstrate  the  effect  on  natural  frequency  of  a  range  of  typical 
soil  stiffnesses.  The  Increase  in  the  first  natural  circular  frequency 
with  soil  stiffness  is  graphed  in  Figure  6.  The  shape  of  the  graph  for 
evidences  the  square-root  relationship  disctissed  earlier  between 
natural  circular  frequency  and  soil  modulus,  E^.  Additional  curves  for 
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the  second  and  third  natural  circular  frequencies  are  merely  sketched  in 
the  figure  to  emphasize  the  diminishing  effect  of  the  soil  stiffness  on 
frequency  for  higher  modes. 

From  these  results,  the  value  f^  =  1.04  hertz  would,  therefore, 

appear  to  underestimate  the  value  corrected  for  the  presence  of  the  sqll 
foundation.  That  is,  the  actual  first  nattiral  frequency  including  the 
effect  of  the  soil  stiffness  could  be  substantially  higher,  perhaps  closer 
to  '2, hertz  according  to  the  graph.  But  these  results  account  only  for 
the  added  influence  of  the  soil  foundation's  stiffness,  and  not  for  soil 
foundation's  inertia.  The  inertial  properties  of  the  soil  cannot  safely 
be  neglected  particularly  for  embedded  structures  such  as  drydocks.  The 
actual  frequencies  may  well  be  similar  to  those  calculated  in  the  absence 
of  soil,  because  the  effect  of  neglecting  soil  inertia  would  appear  to 
mitigate  the  effect  of  neglecting  soil  stiffness.  However,  more  powerful 
soil-structure  interaction  numerical  analysis  models  are  generally  required 
to  also  account  for  soil  foundation  inertia. 

The  preceding  discussion  serves  to  introduce  the  larger  problem 
which  is  the  understanding  of  the  phenomenon  of  soil-structure  interaction. 
Correctly  computing  soil-structure  interaction  effects  when  determining 
the  natural  vibration  and  dynamic  response  of  drydocks  subjected  to  earth¬ 
quake  loads  is  essential.  Both  frequency  domain  and  time  domain  approaches,' 
as  described  by  Wolf  (1985,)  and  Bayo  (198.3,  1987),  need  to  be  considered 
in  future  earthquake  analyses  of  Navy  drydocks  and  waterfront  facilities. 
Linear  and  nonlinear  finite  element  models  of  these  systems  are  generally 
recommended  for  such  analyses. 


TWO-DIMENSIONAL  FINITE  ELEMENT  MODEI^S  OF  DRYDOCKS 

The  purpose  of  this  section  is  to  describe  two-dimensional  finite 
element  structural  analysis  models  of  drydocks  with  particular  emphasis 
on  the  engineering  assumptions  that  are  implied  in  these  models. 

Some  new  calculations  of  the  natural  frequencies  and  natural  mode 
shapes  for  a  two-dimensional  finite  element  model  of  a  drydock  are 
presented  jand  discussed.  These  data  were  prepared  for  a  subsequent 
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compflrison  with  three-dimensional  finite  element  models  which  are  pre¬ 
sented  In  a  later  section.  Two-dimensional  finite  element  models  are 
currently  used  to  support  the  Navy's  drydock  certification  program  with 
respect  to  structural  safety  of  the  drydock.  As  defined  by  Yachnls  (1985) 
drydock  certification  is  an  independent  technical  review  of  existing 
facilities,  and  a  method  for  developing  ways  to  insure  the  safety  of  the 
ship  during  docking  and  while  the  ship  is  being  built  or  repaired.  There 
are  two  points  in  this  definition  that  the  present  discussion  emphasizes: 
(1)  the  procedure  is  fundamentally  a  technical  process,  and  (2)  the  pro¬ 
cedure  encourages  the  Improvement  of,  the  engineering  analysis  methods 
used,  since  the  drydock  certification  program  is  a  continuous  process. 

It  Is'  noted  that  while  the  Naval  Facilities  F.ngineerlng  Command 
(NAVFAC)  has  the  responsibility  for  conducting  certification  studies  of 
drydock  facilities,  the  Naval  Sea  Systems  Command  (NAVSEA)  has  the  ulti¬ 
mate  certification  authority.  Thus,  NAVFAC  is  a  consultant  to  NAVSEA  in 
certifying  the  safety  of  ships  in  drydock.  F>irther,  certification  is  to 
be  conducted  in  accordance  with  MIL-STD- 1625 ,  "Drydocklng  Facilities 
Safety  Certification  Criteria  for  Docking  U.S.  Navy  Ships"  (Naval  Sea 
Systems  Command,  1984). 

According  to  Wu  (1985),  there  are  over  50  graving  drydocks  that 
need  to  be  certified  in  the  United  States.  Naval  shipyards  are  usually 
located  in  areas  with  unfavorable  soil  conditions.  There  have  been  many 
hydraulic  fills  that  utilize  uniform,  fine  sands  which  are  subject  to 
liquefaction.  These  drydocks  were  built  during  the  past  20  to  50  years. 

Drydock  Loads  and  Failure  Modes 

Various  load'  cases  have,  been  defined  for  structural  analysis  of 
drydocks.*  However,  the  earthquake  load  case  for  drydocks  located  in 
seismic  zones  3  and  4  is  of  primary  concern,.  The  earthquake  load  case 
Includes  the  effects  of  dead  load;  ship  blocking  loads;  soil  and  founda¬ 
tion  loads;  and  earthquake- induced,  equivalent  static  loads. 


*It  is  noted  that  blast  loads  have  not  been  ad<lressed. 


General  failure  modes  envisioned  by  Yachnis  (1985)  and  Wu  (1985) 
include:  foundation  instability  where  the  dock  rotates,  shifts,  or 

uplifts  and  generally  displaces  as  a  rigid  body;  structural  damage  to 
the  dock  such  as  overturning  walls;  displacement  of  the  caisson;  failure 
of  utility  tunnels;  and  auxiliary  ifailure  modes  such  as  ship  blocking 
failures  and  overturning  of  cranes  into  the  dock. 

Failure  modes  may  well  be  site-specific  and  include  other  modes  as 
well.  For  example,  burled  pipelines,  water  tunnels,  and  pressure  equali¬ 
zation  systems  which  are  critical  to  operation  and  safety  are  subject  to 
strains  due  to  both  rigid  body  and  relative  motion  (differential  displace¬ 
ment)  of  the  structural  system  (ASCF.  Committee  on  Seismic  Analysis,  1983). 
However,  there  is  evidence  thr.‘.  heavy  industrial  facilities,  provided 
they  have  been  designed  for  seismic  'lOads,  will  perform  very  well  (ASCE 
Committee, on  Dynamic  Analysis,  1987). 

Linear  Static  Analysis  of  Prydooks  for  Earthquake  Toads 

According  to  Chelapatl  and  Takahashl  (1982),  the  earliest  structural 
analyses  conducted  in  support  of  drydock  certification  were  quite  simpli¬ 
fied  and  based  upon  strength  of  materials  theory  (URS/John  A.  Blume  & 
Associates,  1978;  Woodward -Clyde  Consultants,  1979;  and  Moffat  and  Nichol 
Engineers,  1981).  The  general  conclusion  was  that  much  more  can  be  done 
in  applying  available  strtictural  analysis  technology.  Most  of  the  earlier 
studies  did  not  vary  parameters  to  obtain  ranges  of  response,  and  the 
many  assumptions  made  most  likely  lead  to  results  that  are  too  conservative 
It  was  recognized  that  even  a  very  simple,  linear  two-dimensional  finite 
element  model  is  a  more  accurate  basis  for  structural  analysis  than  ele¬ 
mentary  beam  formulas  from  the  theory  of  strength  of  materials. 

Drydock  certification  has  since  been  supported  by  structural  analysis 
methods  which  are  based  on  tvjo-dimens ional  finite  element  modeling  tech¬ 
nology  (Holland  and  Takahashl,  1984;  and  Wu,  et  al.,  1984).  This  method¬ 
ology  for  drydock  analysis  is  reviewed  next. 

Two-dimensional  finite  element  models  of  drydocks  are  based  on  the 
plane  strain  deformation  assumptions  From  the  theory  of  elasticity.  A 
two-dimensional  idealization  of  a  drydock  is  shown  in  Figure  7.  A  typical 
Idealization  of  an  equivalent  static  earthquake  load  is  included  in  the 
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flgiir*.  Tli«*  surrounding  bnr.kfill  sufl  foundstlon  structures  are  excluded 
from  the  model.  Tliese  geotechnical  structures  are  replaced  with  their 
idealized  effects  as  Indicated.  For  example,  the  triangular  distribu¬ 
tions  of  soil  pressure  loads  shown  represent  the  inertia  effect  of' 
adjacent  soil  and  wall  mass. 

Active,  lateral  soil  pressure  is  a.'.siimed  to  increase  linearly  with 
'  depth  according  to  Rankine's  theory  for  the  state  of  s::ress  in  granular 
soil.  The  intensity  depends  on  the  angle  of  friction  and  unit  weight  of 
the  backfill  material  which  varies  according  to  whether  it  is  dry, 
saturated,  or  submerged  (Bowles,  1977).  Ftirther,  a  factor  accounting 
for  wall  movement  may  also  be  introduced  to  increase  the  Intensity. 

There  is  also  a  vertical  component  of  soil  loading  which  is  termed  the 
surcharge,  and  is  repres.inted  in  Figure  ^  by  the  concnni:rated  force 
acting  downward  on  the  sloping  face,  of  the  dr'ydock  wall. 

An  equivalent  static  earthquake-induced  lateral,  soil  pressu’-e  is 
also  Included  in  Figure  7  which  Is  assumed  to  decrease  linearly  with 
depth.  It  is  calculated  by  the  method  discu.ssed  by  Seed  and  Whitman 
(1970).  Tliis  method  is  a  very  simple  extension  of  the  calculation 
method  for  the  static,  active  soil  pressure  mentioned  above.  The  total 
force  (the  area  of  the  active  pressure  triangle)  is  multiplied  by  3/A  of 
the  horizontal  ground  acceleration.  This  product  is  the  Increase  In 
total  lateral  force  due  to  an  earthquake'.  It  is  then  distributed  over 
the  vertical  wall  as  a  trlang\i1ar  pressure  distribution  with  the  centroid 
located  at  2/3  the  wall  height  above  the  drydock  base.  An  additional 
vertical  component  may  be  included'  to  account  for  wall  friction  forces, 
and  where  combined  with  the  lateral  component  the  resultant  equivalent 
static  earthquake- induced  pressure  In  the  drydock  wall  is  obtained  (Wu, 
et  al. ,  1984).  Thus,  on  one  side  of  the  drydock,  a  static  active  soil 
pressure  plus  an  equivalent  static  earthquake-induced  pressure  Is  im¬ 
posed  when  simulating  the  earthquake  load  case. 

On  the  opposite  side,  a  passive  soil  pressure,  which  is  generally 
understood  to  be  greater  than  the  at-rest  active  pressure,  is  imposed. 
Additionally,  a  continuous  column  of  horizontal  springs  is  included  to 
simulate  soil  stiffness  under  earthquake  loading.  The.  statiress  modulus 
is  assumed  to  Increase  linearly  with  depth,  as  shown. 
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the  drydock  floor,  the  hr-ftrlng  stress  of  the  soil  I/’  replaced 
by  a  Winkle  r  foundation  model  which  Is  a  continuous  layer  of  uniformly 
spaced  vertical  springs  with  modulus  K^.  An  elastic  subgrade  modulus  Is 
assigned  to  the  springs  and  Is  determined  according  to  soil  bearing  capacity 
under  assumed  static  load  conditions 

Soil  friction  forces  on  tlie  dryf lock  wall  cannot  bo  calculated  In 
the  current  analysis  procedure  hecause  they  aro  statically  liidetermlnant. 
Various  load  cases  may  he  studied  to  bracket  the  response  and,  thus,  to 
overcome  this  difficulty  Sfwiewhat.  However,  the  basic  problem  of, Indetef 
minacy  remains  because  the  soil  backfill  structure  is  excluded  from  the 
drydock  model;  that  Is,  the  basic  problem  of  so i 1 -structure  Interaction 
Is  simplified  or  decotipled  In  the  present  two-dimensional  stryctural 
analysis  procedure. 

Discussion  of  Hodel  ing  Assiiaipt  ions 

The  current  two-dimensional  finite  element  analysis  procedure,  sup-  , 
porting  drydock  certification  Includes  several  fundamental  assumptions 
as  follows; 

1.  P'lane  strain  ronditlotis  from  elasticity  theory  apply. 

2.  Sol  1 -structure  Interaction  may  he  decoupled. 

.3.  Soil  constitutive  properti.'s  are  linear. 

4.  The  structural  response  Is  independent  of  time. 

These  four  assumptions  are  expanded  and  discussed  below  in  the  order 
listed. 

Plane  Strain  Condition,  Two-dimensional  plane  strain  finite 
ele.ment  analyses  are  often  conducted  for  long,  uniform,  or  prismatic 
striictural/geotechnical  sy.stems  such  as  lined  tunnels  and  earthen  dams. 

A  typical  length-to-width  ratio  of  drydocks  is  .ipprox imately  eight,  which 
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RiAAnN  th*-***  RtnictiirPR  nr#  *iif  f  ir.  ly  long  rtn<l  *lonrl«»r  ««  to  certainly 
anggont  that  plan**  strain  r.ondftlons  tftiilcl  h«»  appropriate.  However,  It 
la  iiaeftil  to  analyze  this  aasiimpt  Ion  iiK)re  t;lo«ely. 

The  conriitloiia  under  whit.h  the  plane  strain  theory  assumption  of 
elasticity  would  be  valid  for  dtydock  analysis  are  stated  here  with  the 
help  o  f  Figure  «: 

a.  The  geometry,  the  loaditig,  and  the  boundary  (support) 
conditions  of  the  drydor.k  dfi  not,  vary  with  coordinate  7.  along 
the  drydork  length. 

b.  Tlte  drydork  displacement  in  the  z-direction,  w,  is  zero 
(i.e.,  the  drydork  does  not  develop  aaial  strain). 

Tliese  conditions,  in  the  context  of  general  two-dimensional,  plane  .strain 
finite  element  analysis  are  discussed  i>y  Callagher  (197'»). 

The  geometry  of  d.ydock  sections  at  a  pumphouse  location  and  at 
stations  near  the  cais.son  and  endwall,  for  example,, are  indeed  different 
from  one  another.  So  the  geometry  of  the  drydock  does  vary  along  its 
length.  Therefore,  the  structural  rigidity  of  a  drydock  sidewall  will 
vary  with  coordinate  7..  Near  the  endwall,  for  instance,  the  huge  In- 
plane  or  membrane  stiffness  of  the  endwall  will  have  a  large  restraining 
effect  on  lateral  displacement  of  the  drydock  sidewall.  Clearly,  results 
from  the  two-dimensional  plane  strain  finite  el'-ment  analysis  will  not 
be  applicable  near  the  endwall.  For  similar  reasons,  the  same  is  true 
near  the  caisson  or  pumphoiise,  or  anywhere  the  drydock  cross-sectional 
geometry  changes  dramatically. 

The  depth  of  the  water  table  aiul  backfill  material  in  many  cases 
will  also  vary  along  the  drydork  lengtii.  Tlierefore,  the  associated  loads 
due  to  hydrostatic  pressure  and  earl  lupiake- ln<lu>-ed  soil  pressures  do 
vary  with  coordinate  7..  It  is  clear  that  condition  (a),  which  is  required 
for  valid  application  of  plane  strain  finite  element  analysis  to  drydocks, 
does  not  strictly  apply.  The  error  which  is  introduced  because  of  this 
must  be  shown  to  be  acceptable.  One  way  to  do  this  would  be  to  conduct 
a  three-dimensional  analysis  and  asse.ss  the  error. 
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Figure  8.  Analysis  of  Plane  Strain  Theory  Assumptions  for  a 


Morpover,  tw6-dlmRns,ionfl  1  flnit»*  olpmcnt  modols  cannot  account  for 
th«  dynamic  torsional  responso  of  structurns  and  for  the  response  for 
nonvert leal 1y  incident  seismic  waves.  The  intrinsic  geometrical  symmetry 
of  two-dimensional  models  raises  serious  doubts  ns  to  their  applicability 
to  three-dimensional  configurations  In  general,  ns  stated  by  Luco  (1982). 
The  approach  to  calculating  a  drydock's  dynamic  response  to  strong-motion 
earthquakes  should,  therefore,  be  based  on  a  three-dimensional  analysis 
If  torsional  and  longitudinal  bending  stresses  are  to  be  considered. 

Soil -Structure  Interaction.  Decoupling  of  the  soi 1 -structure 
interaction  effects  assumes  fundamentally  that  the  structure's  motion  at 
the  soil  interface  does  not  Influence  the  motion  of  the  adjacent  backfill 
and  foundation  materials,  and  vice  versa.  For  example,  this  assumption 
would  neglect  the  affect  of  soil  liquefaction  o,i  the  structure,  and  the 
effect  of  the  structure's  presence  oti  soil  liquefaction.  A. more  accurate 
approach  would  he  to  model  both  the  drydock  and  soil  together  employing 
coupled  soi 1 -structure  interaction  analysis.  This  approach  would  elimi¬ 
nate  the  current  necessity  for  the  many  assumptions  relating  to  soil 
pressure  and  friction  acting  on  the  drydock  walls  and  floor.  Dhatta,  et 
al.  (1985)  show  that  the  modulus  of  the  elastic  springs  that  simulate 
soil  stiffness  next  to  the  structure  under  earthquake  loads  are  a  func¬ 
tion  of  excitation  frequency.  This  frequency  dependence  is  clearly  not 
Included  in  the  present  methodology.  The  properries  of  the  equivalent 
elastic  foiindation  are  assumed  constant  and  are.  based  on  static  soil 
properties.  To  determine  the  significance  of  the  error,  parametric 
studies  varying  the  elastic  spring  moduli  aie  needed.  Alternatively, 
dynamic  analyses  in  the  frequency  domain  may  be  appropriate  (Wolf,  1985) 
where  nonlinear  effect.^  are  minor. 

Soil  Constitutive  Properties.  Invoking  the  assumption  that  founda¬ 
tion  materials  behave  linearly,  forecloses  on  model  predictions  of  failure 
states  within  the  backfill  and  foundation  materials,  or  permanent  plastic 
deformation  of  the  foundation  materials.  Foundation  instability  cannot 
be  predicted  directly  by  the  current  linear  procedure.  To  retain  the 


capability  of  predicting  foundation  failure,  modern  nonlinear  soil  con¬ 
stitutive  models  must  be  used  for  the  backfill  and  foundation  structures. 
The  Importance  of  a  properly  operating  backfill  relief  system  is  empha¬ 
sized  by  Zola  and  Boothe  (1960)  for  fully  re]leve<l  drydock  designs. 

Structural  Response.  The  assumption  that  the  structure's  response 
is  independent  of  time  ignores  the  possibility  of  resonance  or  amplified 
response  in  the  drydock  when  the  earthquake  excitation  frequency  matches 
the  drydock  system's  natural  frequency.  Resonance  effects  would  seem  to 
be  particularly  Important  to  the  dynamic  response  of  blocked  ships,  dry- 
dock  cranes,  and  other  similar  "superstructures"  in  contact  with  the 
drydock  itself.  These  considerations  bear  directly  on  ship  safety, 
which  is  the  primary  goal  of  the  drydock  certification  program. 

Other  Considerations.  Other  consideration.s  which  foretell  the  need 
for  dynamic  analysis  of  drydocks  are  as  follow.s.  Amplification  of  ship 
blocking  loads  could  be  calculated  by  a  dynamic  analysis  procedure.  It 
is  neglected  in  the  current  procedure.  A  dynamic  analysis  would  capture 
the  effects  of  radiation  damping  in  the  soil  foundation,  which  is  cur¬ 
rently  bellaved  to  be  an  important  energy  dissipation  mechanism  during 
dynamic  events.  This  reinforces  the  need  for  dynamic  soil-structure 
interaction  models. 

The  technical  problems  associated  with  current  two-dimensional  , 
finite  element  models  of  drydocks  are  summarized  in  Table  1.  A  correct 
procedure  for  earthquake  analysis  of  drydocks  would  be  based  upon  a  three 
dimensional,  nonlinear  dynamic  analysis  of  the  drydock  including  the 
surrounding  backfill  and  foundation  structures. 

Natural  Frequencies  and  Mode  Shapes 

Future  dynamic  analysis  of  drydocks  would  necessarily  be  preceded 
by  an  eigenvalue  analysis  of  the  finite  element  model  to  determine  the 
natural  frequencies  and  the  natural  mode  shapes  of  the  drydock.  This 
section  presents  and  discusses  the  results  of  such  an  analysis  for  a 
two-dimensional  model.  The  two-dimensional  finite  element  model  in  this 
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Tablri  1.  Technical  Problems  With  Current  Two-Dimensional 
Finite  Element  Analysis  of  Drydocks 


1.  2-D  PLANE  STRAIN  CONDITIONS  DO  NOT  APPLY 

•  Torsional  and  longitudinal  bending  response 
cannot  be  replicated. 

2.  SOIL-STRUCTURE  INTERACTION  IS  NOT  INCLUDED 

•  Drydock  and  foundation  response  is  inaccurate. 
SOIL.  CONSTITUTIVE  PROPERTIES  ARE  NOT  LINEAR 

•  Foundation  material  failure  cannot  be  detected. 
4.  STATIC  ANALYSIS  CANNOT  DETECT  RESONANCE 

•  Predicted  response  may  not  be  conservative. 


case  is  a  one-element-thick  slice  of  a  three-dimensional  model  of  dry- 
dock  No.  6,  Puget  Sound  NSY.  Natural  frequencies  and  mode  shapes  for  a 
three-dimensional  model  of  this  drydock  will  also  be  presented  and  dis¬ 
cussed  in  a  subsequent  section  of  this  report.  Results  for  the  two'- 
dimensional  model  and  the  three-dimensional  model  will  be  compared  to 
determine  differences  in  the  fiindamental  dynamic  properties  of  the  two 
modeling  approaches. 

To  obtain  a  two-dimensional  model  with  plane  strain  behavior,  dis¬ 
placements  in  the  long  direction  of  the  drydock  for  the  one-element-thick 
three-dimensional  model  are  constrained  to  be  r.p.ro.  This  imposes  con¬ 
dition  (b)  for  plane  strain  behavior  given  in  the  previous  section. 
Condition  (a)  is  tacitly  assumed  to  hold.  Thus,  the  model  presented 
here  and  used  for  calculating  frequencies  and  mode  shapes  would  be 
equivalent  to  a  natural  vibration  analysis  of  drydock  No.  6  based  on  the 
current  two-dimensional  plane  strain  modeling  practice. 

In  this  case,  the  two-dimensional  finite  element  model  exploits 
symmetry  about  the  centerline  of  the  drydock,  as  shown  in  Figure  9.  A 
consequence  of  this  is  that  all  natural  frequencies  and  mode  shapes 
associated  with  modes  of  vibration  that . are  ant isymmetrlcal  about  the 
centerline  will  not  be  admitted  in  the  calculation.  Only  natural  mode 
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shapes  that  are  symmetrical  with  respect  to  the  drydock  centerline  will 
result  from  this  analysis.  To  obtain  all  modes,  symmetrical  and  anti- 
symmtetrical,  would  be  straightforward,  but  a  full  two-dimensional  slice 
of  the  drydock  would  have  to  be  modeled.  Nonetheless,  the  results  from 
the  present  symmetrical  two-dimensional  model  are  a  sufficient  basis  for 
subsequent  comparison  with  the  results  from  both  one-dimensional  and 
three-dimensional  models. 

Other  features  of  the  model  are  indicated  in  Figure  9(b).  A  single 
row  of  soil  elements  represent  the  effect  of  the  surrounding  soil  founda¬ 
tion.  The  drydock  material  model  and  foundation  material  model  are  assumed 
to  be  Isotropic  and  linear  elastic.  An  equivalent  modulus  of  elasticity 
and  Poisson  ratio  of  the  foundation  were  calculated  from  a  typical  value 
for  soil  bearing  modulus. 

The  first  five  symmetrical  natural  mode  shapes  and  frequencies  are 
shown  in  Figures  10(a)  through  10(6).  Only  the  concrete  drydock  portion 
of  the  model  is  shown  in  these  figures;  the  utility  tunnel  and  the  row 
of  soil  foundation  elements  are  omitted  from  the  graphics  presented. 

These  results  show  that  the  first  mode  shape  is  a  floor  bending 
mode  and  the  second  mode  shape  is  a  cantilevered  wall  bending  mode. 
Subsequent  modes  are  higher  order  combinations  of  these  two  deformation 
mode  shapes.  The  equivalent  static  analysis  procedure  discussed  in  the 
previous  section  effectively  assumes  a  deformation  mode  shape  that  is 
similar  to  a  combination  of  these  first  two  mode  shapes.  Thus,  it  may 
be  concluded  that  the  current  equivalent  static  analysis  procedure 
effectively  neglects  any  contribution  from  modes  of  deformation  beyond 
the  second. 

The  lowest  natural  frequency  calculated  for  the  two-dimensional 
plane  strain  finite  element  model  is  2  hertz.  This  might  appear  high 
compared  to  that  which  was  calculated  for  the  one-dimensional  models  in 
the  previous  section,  which  was  abotit  1  hertz.  However,  the  values  are 
not  directly  comparable  because  they  are  fundamentally  different  mode 
shapes  from  fundamentally  different  models.  In  the  present  case  they 
are  lateral  vibration  modes,  while  in  the  former  case  they  are  longitu¬ 
dinal  vibration  modes.  In  a  dynamic  analysis  with  an  input  ground  motion 
having  a  prominent  1-hertz  energy  component,  use  of  this  two-dimensional 
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model  would  lead  to  results  having  minimal  resonance  behavior.  However, 
the  one-dimensional  model  would  respond  with  strong  resonance.  Thus, 
depending  on  which  model  was  being  used,  dramatically  different  conclu¬ 
sions  on  the  Importance  of  basic  dynamic  behavior  could  be  drawn  from 
their  results.  The  one-dimensional  model  considers  only  longitudinal 
bending  modes  and  the  two-dimensional  model  considers  only  wall  and 
floor  bending  modes  in  a  transverse  plane. 

Another  observation  from  these  results  is  that  the  wall  bending 
modes  are  high-frequency  modes  (l.e.,  greater  than  2  hertz).  This  is  a 
result  of  the  very  large  bending  stiffness  of  the  wall  and  floor  sections, 

particularly  where  they  meet  at  the  base  of  the  vertical  cantilever.  As 

.  / 

a  result,  the  walls  would  rarely  be  expected  to  bend  in  modes  higher 
than,  for  example,  the  third  because  most  strong  motion  earthquakes  do 
not  contain  strong  energy  components  beyond  5  hertz.  The  first  two  modes 
of  wall  deformation  are  effectively  subsumed  by  the  present  equivalent 
static  analysis  procedure  for  earthquake  loads.  Thus,  the  equivalent' 
static  procedure  may  be  a  good  approximation-  to  the  two-dimensional 
dynamic  analysis.  However,  the  sidewalls  will  also  bend  in  longitudinal 
vibration  modes  which  are  of  course  precluded  in  the  two-dimensional 
model  approach. 

TVo-dimensional  finite-element  models  of  drydocks  normally  include 
a  Winkler  model  (a  row  of  linear  elastic  springs)  for  representing  the 
effect  of  the  soil  foundation.  To  study  the.  difference  between  the 
effects  of  using  a  row  of  two-dimensional  finite  elements  and  a  row  of 
one-dimensional  springs,  the  latter  model  was  also  considered  in  this 
study.  The  first  ten  natural  frequencies  and  mode  shapes  using  an 
equivalent  Winkler  model  for  the  drydock  foundation  were  calculated. 

The  results  for  both  soil  foundation  models  were  very  similar.  The  mode 
shapes  were  basically  the  same  as  those  presente.d  in  Figure  10  and, 
therefore,  are  not  presented.  The  natural  frequencies  were  also 
similar,  particularly  for  the  more  important  lower  modes.  The  natural 
frequencies  for  the  first  ten  modes  calculated  with  both  foundation 
models  are  graphed  in  Figure  11.  Tlius,  it  does  not  seem  to  matter 
greatly  which  type  of  model  (Winkler  or  finite  element)  is  used  tq 
account  for  the  foundation  when  considering  the  use  of  two-dimensional 
finite  element  models  of  drydocks. 
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TWO-DIMENSIONAI.  LUMPED  PARAMETER  MODEL  OF  BI/50KED  VESSEL 


A  principal  goal  of  the  Navy  Drydock  Certification  Program  is  the 
establishment  of  the  safety  of  the  docked  vessel.  One  important  aspect 
that  can  arise  is  consideration  of  the  safety  of  a  blocked  vessel 
relative  to  earthquake  loading. 

The  need  for  earthquake  analysis  of  drydocks  has  been  questioned 
technically  on  the  basis  that  a  drydock  is  an  embedded  structure  having 
no  superstructure,  as  compared  with  a  high-rise  building,  for  example. 
However,  a  blocked  vessel  may  be  considered  a  subsystem  of  the  entire 
drydock  system,  no  differently  than,  for  example,  a  pumphouse  subsystem 
or  a  caisson  subsystem.  In  this  sense,  a  blocked  vessel  may  Indeed  be 
viewed  as  a  "superstructure"  that  responds  dynamically  to  input  base 
motion  from  the  drydock  floor.  Further,  the  analysis  of  a  complete 
drydock  system,  blocked  vessel  included,  is  in  principle  within  the 
capability  of  modern  structural  dynami  analysis  technology.  There  are 
various  modeling  approaches  that  could  be  applied  in  this  regard,  some 
more  complex  and  possibly  more  accurate  than  others.  However,  it  is 
believed  that  the  dynamic  interaction  of  a  blocked  vessel  with  the 
drydock  floor  could  be  calculated  with  mod  rn  methods  to  obtain  accurate 
data  on  the  dynamics  of  the  vessel,  blocking  system,  and  drydock  floor. 
The  purpose  of  this  section  is  to  demonstrate  this  using  a  simple,  two- 
dimensional  lumped  parameter  model. 

Construction  of  a  fully  three-dimensional  subsystem  model  of  a 
blocked  vessel  is  possible.  However,  for  the  present  demonstration,  a 
two-dimensional,  lumped  parameter  model  of  a  blocked  submarine  was 
developed  and  analyzed  for  its  dynamic  response  relative  to  a  simulated 
drydock  floor  excitation.  The  model  used  here  was  inspired  by  Hepburn, 
et  al.  (19a8),  a  study  which  was  sponsored  by  NAVSEA.  This  publication 
also  gives  an  account  of  the  concern  for  the  safety  of  blocked  Navy  . 
vessels  relative  to  earthquake  loads.  Tt  develops  a  considerable  amount 
of  engineering  data  demonstrating  that  a  blocked  vessel  model  can  be 
analyzed  using  the  methods  of  structural  dynamics  to  obtain  accurate 
blocking  system  design  data.  The  importance  of  knowing  aprlorl  the 
floor  excitation  Imposed  by  the  drydock  system  is  also  emphasized.  To 
the  extent  that  a  drydock  responds  to  an  earthquake  as  a  deformable  body 
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rather  than  as  a  rigid  body,  it  will  modify  the  ground  motion  to  which 
the  blocked  vessel  is  otherwise  exposed.  Specification  of  this  modified 
ground  motion  is  required  for  successful  dynamic  analysis  of  the  blocked 
vessel.  This  is  achievable  by  conducting  soil-structure  interaction 
analyses,  as  discussed  earlier  in  this  report. 

Lumped  Parameter  Model  of  a  Blocked  Submarine 


The  drydock  is  assumed  to  be  a  rigid  body  that  excites  the  blocked 

vessel  simultaneously  in  both  a  horizontal  and  vertical  direction.  The 

lumped  parameter  model  is  shown  in  Figure  12  where  a  cylindrical  mass, 

representing  a  submarine,  is  spring-mounted  to  the  drydock  floor  and 

walls.  In  a  lumped  parameter  model  the  system's  inertia  and  elastic 

properties  are  assumed  to  be  concentrated  at  discrete  points.  For 

simplicity,  the  submarine  is  assumed  to  be  a  uniform  solid,  circular 

cylinder.  It  has  radius,  r,  mass,  m,  and  mass  moment  of  Inertia  about 

its  center  of  gravity,  I  .  The  blocking  system  Is  assumed  to  consist 

eg 

of  keel  blocks,  side  blocks,  and  horizontal  braces  representing  a  pro¬ 
posed  wale  shoring  system  consisting  of  rubber-tipped  steel  beams  canti¬ 
levered  out  from  the  drydock  sidewalls.  These  support  elements  are 
represented  by  linear  elastic  springs  having  moduli  kj^  through  k^,  as 
shown. 

In  this  case,  the  lumped  parameter  model  has  three  degrees  of 
freedom.  The  variables  u,  v,  and  0  represent  the  absolute  displacement 
components  of  the  submarine's  center  of  gravity  in  the  horizontal  and 
vertical  directions,  and  the  rotational  displacement  about  the  renter  of 
gravity,  respectively.  Instead  of  absolute  displacement  components,  u 
and  v,  the  displacements  of  the  center  of  gravity  relative  to  the  dry- 
dock  may  be  used.  These  are  indicated  by  the  variables  x  and  y  (see 
Figure  12).  The  relative  displacement  components  are  generally  more 
useful  when  t.he  excitation  is  base  motion.  Therefore,  the  d  ispl aoemen* 
variables  used  In  this  demonstration  are  x,  y,  and  ^  Tie  prescribed 
Input  base  motion  Is  represented  by  Xj^  and  y^^  wljj'r  4'*-  essrwied  to  be 
simple  harmonic  functions  representing,  respective;/,  th»  bc-rizonta'!  and 
vertical  displacements  of  the  drydock. 


34 


Rgure  12.  Two-Otmensional  Lumped  Parameter  Model  of  a  Blocked  Submarine  in  Orydock. 

The  values  used  for  the  submarine’s  ass  and  rotational  inertia  are 
calculated  based  upon  a  displacement  of  3,570  long  tons.  Fourteen  iden¬ 
tical  blocking  stations,  of  the  type  described,  were  assumed  to  support 
the  vessel  along  its  length.  Tims,  the  values  iised  for  m  and  I  repre- 

»  I 

s<-nT.  physical  properties  "per  blocking  station,"  and  are  suitable  for 
use ' in  the  proposed  two-H imens 'ona 1  or  planar  model  (Figure  12). 

The  elastic  properties  ii.sed  for  the  blocking  system  were  calculated, 
In  some  cases  roughly,  from  data  presented  bv  Hepburn,  et  al.  (19B8). 
Individual  blocks  are  often  composed  of  mn I i iplo  materials  including 
concrete,  rubber,  and  ork  and  fir  timber  components.,  The  methods  for 
cal  dilating  equivalent  elastic  spring  mod’nli,  especially  for  the  lateral 
sprtftg  «oduli  of  the  side  and  keel  blocks,  can  be  complex  when  accounting 
"oth  sheer  and  axial  de form.at ion ,  as  well  as  for  the  many  materials 
; -..ad ,  Yet,  accurate  analytical  results,  fio  depend  on  accurate  input 
dfSP'  these  ami  other  elastic  properties  of  the  model.  The  use  of 

'epsr  i'o'ator  blocking  systems  (i.c.,  elastomeric  bearings)  would 

•wse  c  '  S' e  this  uncertainty,  hut  otherwi -i  field  measurement  program 


Id  needl'd  to  obtain  rnasonably  acciirnto  data  for  tho  nlastiC'  proportlas 
of  ciirrpnt  blocking  materiald  and  sydtpms.  ValiiPd  u»pd  for  the  parameters 
of  the  present  lumped  parameter  model  are  listed  in  Table  2,  and  are 
sufficiently  accurate  for  this  demonstration. 


Table  2..  Parameters  for  Blocked  Submarine  Model 


*^1 

=  5.30  X  10^  Ib/ln. 

r  = 

1«0  in. 

=  5. 10  X  10*  Ib/in. 

Id  = 

i.4«0  X  10^  In-s^/in. 

=  3.0  X  10^,  Ib/in. 

I  / 
rfi 

=  2.397  X  10^  Ib-s^-ln. 

=  6.0  X  10*  Ib/in. 

a  - 

30** 

Equations  of  Motion 

The  equations  of  motion, for  the  lumped' parameter  model  are  con¬ 
veniently  derived  by  applying  either  Newton's  Mws  or  Lagrange's  equa¬ 
tions  of  motion  to  the  model  of  Figure  12.  The  results  are  a  set  of 
three,  second  order  differential  equations  in  time.  In  matrix  form,  the 
equations  of  motion  are; 


Dt 

0 

0 

X 

2(k^ 

4k3) 

0 

-2k3C 

X 

-mx, 

b 

0 

m 

0 

y 

+ 

0 

(kj  +  2k2) 

0 

y 

= 

-myb 

0 

0 

I 

0 

-2k,C 

0 

2fk,C^  +  k.,S^) 

0 

. 

cgj 

.  . 

3 

3  2  ■ 

0 

or, 

m  X  +  k  X  =  3 

The  right-hand  side  term  Q  is  the  force  vector  .and  it  Includes  the  terms 
-mXj^  and  which  are  referred  to  as  effective  forces.  They  represent 

the  system's  excitation  for  any  prescribed  ground  acceleration  components, 
Xj^  and  y^^.  The  terms  C  and  S  denote  the  constants  r  cos  a  and  r  sin  a, 
respectively,  which  define  the  {josition  of  the  side  blocks. 
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In  this  dnmonstratlon,  a  strong-mot  I r>ii  earthquake  is  simulated  by 
specifying  the  ground  accelerations  to  be  a  sln'isotd  as  follows: 

X|^  =  y^  =  0.25g  sln(2iT/T)t 

where  g  is  the  acceleration  of  gravity,  and  T  is  the  period  of  the  simple 
harmonic  motion.  A  value  of  1  second  was  used  for  T.  This  excitation 
simulates  a  substantial  earthquake  with  a  0.25g  peak  acceleration  and  a 
frequency  of  1  hertz. 

In  the  equations  of  motion,  it  is  noted  that  the  second  equation  is 
uncotipled  and  may  be  solved  independently.  This  is  due  to  the  tacit 
assumption  of  small  displacements  in  the  blocked  submarine's  response, 
and  the  symmetry  of  the  blocking  system.  Motion  in  the  y-dlrection  can 
exist  independently  of  motion  in  the  other  directions.  Further,  because 
the  forcing  function  is  elementary,  this  equation  can  be  solved  analyti¬ 
cally.  That  leaves  the  other  two  equations  of  motion  to  be  solved  simul¬ 
taneously,  because  they  are  elastically  coupled.  Motion  of  the  submarine 
in  the  x-dlrection  cannot  exist  without  motion  in  the  0-dlrectlon,  and 
vice  versa,  according  to  this  model.  Despite  the  coupling,  the  remaining 
two  equations  can  also  be  solved  analytically.  For  more  complex  formula¬ 
tions,  the  equations  of  motion  can  he  solved  numerically  using  either 
the  mode  superposition  solution  method  or  the  direct  numerical  integration 
solution  method,  as  described  by  Craig  (1981),  for  example. 

Natural  Frequencies  and  Natural  Mode  Shapes 

Before  determining  the  dynamic  response,  it  is  useful  to  determine 
the  fundamental  natural  frequencies  and  natural  mode  shapes  of  the  model 
to  enhance  understanding  of  its  dynamic  properties  and  behavior.  These 
data  are  obtained  by  solving  the  eigenvalue  problem: 


for  the  natural  circular  frequencies  w  and  the, natural  mode  shapes  (J. 

Since  there  are  three  degrees  of  freedom  (thred  displacements),  the  system 
correspondingly  has  three  frequencies  and  three  mode  shapes.  It's  worth 
noting  that  the  eigenvalue  solution  can  also  be  obtained  analytically 
since  only  two  of  the  three  equations  are  coupled  In  the  above  matrix 
eigenvalue  equation.  Otherwise,  a  microcomputer  system  may  be  used  with 
a  microcomputer  program  such  as  CAI.80  (Bayo  and  Strubb,  1988),  which  has 
the  capability  to  solve  the  eigenvalue  problem  as  well  as  to  solve  for' 
the  dynamic  response  of  small  to  medium  structural  models.  In  this  case, 
both  methods  were  used  to  verify  accuracy  of  the  computer  solution  and 
both  methods  yielded  virtually  identical  solutions. 

The  results  of  the  calculated  natural  frequencies  and  natural  mode 
shapes  are  presented  in  Figure  13.  The.  first  calculated  natural  frequency 
is  1.01  hertx,,and  therefore  a  substantial  resonance  In  the  first  mode 
of  vibration  will  occur  In  the  dynamic  response  since.  In  this  demonstra¬ 
tion,  the  prescribed  excitation  frequency  was  chosen  to  be  1  hertz.  The 
first  mode  shape  sketched  is  coupled,  but  its  vibration  Is  primarily 
rotational.  This  is  a  result  of  the  very  large,  rotational  inertia  of 
the  submarine. >  The  sketch  of  the  mode. shapes  indicates  the  coupling 
that  exists  between  horizontal  vibration  and  rotational  vibration  In 
both  modes  1  and  2,  whereas  mode  3  is  uncoupled  and  consists  of  vibration 
solely  in  the  vertical  direction.  Thus,  it  may  be  expected  that  for  an 
earthquake  which  tends  to  cause  horizontal  floor  motion,  the  blocked 
submarine  will  have  a  tendency  to  vibrate  rotat tonally  about  its  center 
of  gravity  in  addition  to  vibrating  horizontally.  However,  for  an  earth¬ 
quake  that  tends  to  cause  vertical  floor  motion  primarily,  the  blocked 
submarine  wi 11  tend  to  vibrate  primarily  in  the  vertical  direction. 

Dynamic  Response  to  Simplified  F.arthquake  T.oad 

To  obtain  the  dynamic  response  of  the  blocked  submarine  subjected 
to  the  simplified  earthquake,  the  mode  superposition  method  of  solution 
was  employed.  The  solution  was  calculated  using  the  CALSO  computer 
program  and  also  by  hand  as  an  analytical  check  on  the  computer  solution 
of  the  dynamic  response.  The  two  dynamic  response  results  were  again 
virtually  Identical. 
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Natural  Frequencies: 


39 


Rgure  13.  Natural  Frequencies  and  Mode  Shapes  for  a  Blocked  Submarine  in  D.rydock. 


The  dynamic  response  restilts  for  the  three  displacement-time 
histories  x(t),  y(t)  and  0(t)  are  shown  graphed  for  a  duration  of  almost 
4  seconds  in  Figure  14(a).  The  drydock  floor  excitation  is  juxtaposed 
with  the  response  data  in  Figure  14(b)-  The  results  show  the  rotational 
vibration  growing  due  to  a  resonance  condition  in  the  absence  of  any 
assumed  damping  components  in  the  blocking  system  as  modeled.  This  is 
essentially  the  mode  1  resonance  condition  which  was  anticipated,  and  it 
predicts  large  rotational  oscillations.  For  example,  after  shaking  for 
3  seconds  the  rotation  amplitude  has  increased  to  0.02  radians  which, 
amounts  to  about  3.6  Inches  of  tangential  displacement  on  the  submarine 
hull's  surface.  Also,  during  the  first  3  seconds  of  shaking,  the  maxi¬ 
mum  displacements  of  the  hull  center  of  gravity , relatlye  to  the  drydock 
are,  respectively,  3.5  inches  and  2.5  inches,  in  the  horizontal  and 
vertical  directions. 

Though  an  effort  was  made  to  make  the  model  parameters  in  this 
demonstration  as  realistic  as  possible,  these  results  are  meant  only  to 
illustrate  how  tlie  dynamic  characteristics  of  a  blocked  vessel  during 
drydock  shaking  can  be  calculated.  The  interaction  between  the  vessel 
and  the  drydock  is  important  to  the  understanding  of  the  force  and 
displacement  effects  for  the  vessel  as  well  as  the  dynamic  forces  acting 
on  the  drydock  floor;  More  attention  must  be  given  to  an  accurate  des¬ 
cription  of  the  prescribed  model  parameters,  and  the  prescribed  accelera¬ 
tion  of  the  drydock  floor. 

Measured  Earthquake  Response  of  Drydock 

On  1  October  1987,  a  magnitude  5.9  earthquake  with  a  0.45g  maximum 
acceleration  occurred  at  Whittier  Narrows  in  the  bos  Angeles  Basin.  The 
epicenter  was  approximately  42  km  to  the  northeast  from  drydocks  1,  2, 
and  3  sited  at  the  Long  Beach  NSY.  Figure  15  shows  the  relative  loca¬ 
tions  of  these  drydocks  to  each  other.  Strong-motion  accelerograph 
instrumentation  placed  at  two  points  on  both  drydocks  1  and  2  recorded 
the  responses  of  the  drydocks. 
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(a)  Dynamic  Response 
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(b)  Excitation 


Figure  14.  Dynamic  Response  of  a  Biocked  Submarine  in  Drydock, 


The  USS  LEAHY  (CG-16)  was  blocknd  in  drydock  3  during  the  event, 
and  sustained  sufficient  motion  to  cause  permanent  shifting  of  the  hull 
relative  to  the  side  blocks.  The  relative  displacements  were  recorded 
photographically,  and  discussed  In  the  paper  by  Hepburn,  et  al.  (1988). 

It  Is  believed  the  drydock  sustained  only  a  0,05g  maximum  acceleration 
during  the  event.  Hepburn,  et  al.  (1988)  also  describes  a  study  In 
which  post-predictions  from  a  lumped  parameter  model  of  the  blocked 
vessel  successfully  matched  thd  recorded  pattern  of  blocking  movement. 

The  purpose  of  this  section  Is  to  provide  a  summary  analysis  and  dis¬ 
cussion  of  the  corrected  time  histories  of  motion  from  the  accelerograph 
traces  which  were  compiled  at  the  Naval  Civil  Engineering  Laboratory 
(NCEL)  by  Lew  (1988),  and  used  in  the  model  study  by  Hepburn,  et  al. 
(1988). 

Accelerographs  recorded  the  motion  of  drydock  1  at  two  stations, 

580  feet  and  875  feet  from  the  caisson,  along  the  length  of  the  drydock. 

The  time  histories  of  the  transverse  component  of  the  motion  at  the  two 

locations  are  presented  In  Figures  16(a)  and  16(b).  The  peak  values  of 
the  motion  are  Indicated  in  these  figures.  While  the  peak  values  of  the 
motion  at  the  two  locations  are  not  equal,  they  do  occur  at  the  same 
time.  Further,  the  time  histories  for  the  two  locations  are  remarkably 
similar.  This  can  be  observed  by  juxtaposing  the  two  sets  of  graphs  in 
Figure  16(a)  and  16(b).  Two  points  295  feet  distant  had  virtually  the 
same  kinematics  in  the  transverse  direction.  The  similarity  is  strongest 

for  the  acceleration.  Some  deterioration  in  this  similarity  is  noted 

for  the  velocity  and  displacement  (this  is  expected  from  the  smoothing 
Inherent  In  single  and  double  integration  of  corrected  acceleration  data, 
from  which  these  graphs  are  obtained).  These  data  are  clear  evidence 
that  drydock  1  responded  essentially  as  a  rigid  body  rather  than  as  a 
deformable  body  In  the  transverse  direction  during  the  Whittier  Narrows 
earthquake. 

Similar  conclusions  can  be  drawn  from  an  analysis  of  the  data  for 
the  vertical  and  longitudinal  components  of  the  motion  of  drydock  1. 

Thus,  It  Is  believed  that  drydock  1,  Long  Beach  NSY,  responded  as  a 
rigid  body  to  the  1  October  1987  Whittier  Narrows  earthquake. 
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PEAK  VALUES:  ACCELERATION  >  -34.07  CM/SEC/SEC  VELOCITY  »  5.06  CM/SEC  DISPLACEMENT  »  -1 .21  CM 


6.  Whittier  Narrows  Earthquake  -  Induced  Motion  of  Dtydock  1 ,  LBNSY,  Transverse  Direction. 


PEAK  VALUES:  ACCELERATION  =  -41 .62  CM/SEC/SEC  VELOCITY  =  6.05  CM/SEC  DISPLACEMENT  »  -1.43  CM 


Rgure  16.  Continued. 


The  data  recorded  for  drydock  2  were  neither  as  extensive  nor  as 
conclusive.  There  was  some  correlation  between  the  longitudinal  and 
transverse  acceleration  time  histories.  Again,  correlation  diminishes 
for  velocity  and  displacement.  There  was  no  apparent  correlation 
between  the  vertical  acceleration  time  history  and  the  transverse  or 
longitudinal  time  histories. 

An  estimate  of  the  natural  period  of  the  motion  for  drydock  1  can 
be  obtained  by  the. number  of  zero  crossings  In  the  displacement  graphs 
of  Figure  16  for  a  specified  duration.  The  result  Is  an  estimated 
natural  period  of  1  second,  or  an  estimated  natural  frequency  of  1 
hertz,  for  the  excitation  of  drydock  1. 

Since  It  Is  believed  that  this  is  rigid  body  motion,  this  natural 
frequency  must  also  exist  for  the  soil  system  In  which  the  drydock  is 
founded.  At  this  particular  site,  the  soil  is  soft  hydraulic  fill  and 
beneath  this  are  alternating  layers  of  silts,  sands,  and  clays  which  In 
turn  overlay  an  aqnlfet,  i^cale  model  tests  of  the  effects  of  embedment 
length  on  simple  fon-.i-.’at  > :>ns  conducted  by  Jennings  and  Wu  (1984)  tend  to 
show  the  foundation  ar<^i  s  ’il  move  as  a  rigid  body  for  deep  foundations. 

Actual  field  data  1  )r  the  response  of  a  large  Navy  drydock  to  earth¬ 
quake  loads  are  valuaile.  In  this  case,  the  drydock  response  is  believed 
to  be  rigid  body  motio,,  ard  the  frequency  of  motion  is  influenced  by  the 
soil  column  upon  whid  d.- /dock  1  Is  sited.  The  Input  acceleration  is, 
however,  remarkably  snail,  and  does  not  qualify  as . strong-motion  accelera¬ 
tion.  Thus,  the  meaniiig  of  these  data  Is  limited  In  regard  to  the  response 
of  drydocks  to  strong-rtt,;t  ion  earthquakes.  Nonetheless,  however  small 
the  excitation,.  It  was  sufficient  to  cause  damage  to  the  blocking  system 
supporting  a  ship. 

THREE-DIMENSIONAL  FINITE  EHiMENT  MODEI^  OF  DRYDOCKS 

Solids  Modeling 

Construction  of  three-dimensional  finite  element  models  requires 
substantial  computer  resources  almost  equal  to  those  required  for  the 
actual  analysis  of  the  model.  Computing  power  and  computer  storage  are 
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required  to  manipulate  and  maintain  large  data  files  of  information  on 
geometry,  large  numbers  of  nodes,  and  the  connectivity  or  relationship 
of  elements  to  one  another.  As  the  data  base  grows,  simple  commands 
like  inserting  a  single  nodal  point  during  model  generation  may  take 
minutes,  much  to  the  chagrin  of  the  user  who  is  constructing  the  model. 

The  ease  of  constructing  three-dimensional  finite  element  models  is 
dependent  on  the  type  of  geometric  model  needed.  Three-dimensional  geo¬ 
metric  models  fall  into  three  categories:  (1)  wireframes,  (2)  surfaces, 
and  (3)  solids.  A  brief  discussion  on  how  each  pertains  to  finite 
element  modeling  is  presented  in  Appendix  B. 

Substructure  Model lug 

Managing  the  complexity  of  a  large  three-dimensional  finite  element 
analysis  of  a  drydock  system  could  be  greatly  facilitated  by  substructure 
modeling  procedures.* 

Substructure  theory  is  a  generalization  of  the  static  condensation 
procedure  in  matrix  structural  analysis  or  in  the  finite  element  method 
where,  for  example,  the  nodal  variables  for  an  unattached  internal  node 
of  a  quadrilateral  element  are  eliminated  in  terms  of  the  nodal  variables 
for  the  four  corner  nodes.  In  substructure  modeling,  the  selection  of 
the  internal  variables  or  degrees  of  freedom  (dof)  which  are  to  be 
eliminated  is  arbitrary,  so  that  a  large  structiiral/geotechnlcal  system 
may  be  subdivided  into  substructures  arbitrarily,  and  the  dof  internal 
to  each  substructure  may  then  be  eliminated,  while  retaining  only  the 
dof  at  the  substructure  interface  level.  The  theory  of  substructure 
modeling  in  the  context  of  drydock  systems  is  reviewed  in  Appendix  C. 

Substructure  modeling  is  but  one  approach  to  overcoming  limited 
computer  resources  (speed  and  memory)  when  working  with  very  large 
structural  models.  For  effective  use,  analysts  should  be  very  familiar 
with  the  substructure  modeling  source  code.  Substructure  modeling, 
while  simple  in  theory,  is  difficult  to  implement  in  large-scale  com¬ 
puter  programs.  A  certain  degree  of  substructure  capability  is  Included 

*Przemienckl  (1968)  provides  a  classic  introduction  to  the  theory  of 
substructures,  while  Furulke  (1972)  provides  a  more  topical  treatment 
of  the  method  in  the  context  of  large-scale  finite  element  problems. 
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In  some  commercial  computer  programs  such  ns  ADfNA  (AI)INA  Engineering, 
Inc.,  I9R5),  ABAQUS  (HLbbltt,  Karlsson  and  Sorensen,  Inc,,  1988),  and 
GIFTS  (CASA/GTFTS,  Inc.,  1988). 

In  practical  application,  there  are  three  good  reasons  why  sub¬ 
structure  modeling  Is  employed: 

1.  The  same  substructure  repeats  Itself.  When  many  substructures 
are  Identical,  the  method  becomes  computationally  efficient. 

2.  Nonlinear  behavior  confined  to  a  small  part  of  the  structure. 
Computational  efficiency  derives  from  reducing  the  size  of 
expensive  nonlinear  calculations. 

3.  Organization  of  large-scale  modeling  of  complex  structures. 

For  example,  office  A  addresses  the  concrete  drydock,  office  B 
addresse?  the  pump  house,  and  office  C  addresses  the  steel 
caisson. 

Drydock  Substracture  Models 

The  construction  of  a  three-dimensional  finite  element  model  of 
drydock  6  located  at  the  Puget  Sound  NSY  was  attempted  based  upon  the 
Idea  of  substructure  modeling.  The  model  was  to  be  built  using  eight- 
node  solid  brick  elements  for  the  walls  and  floor  of  the  drydock,  four- 
node  plate  elements  for  the  utility  tunnels  atop  the  sidewalls,  and 
elastic  springs  for  the  Jitility  tunnel  supports  and  surrounding  soil. 

The  material  properties  of  the  concrete  and  the  steel  rebar  would  be 
"smeared"  to  yield  an  effective  material  stiffness  for  the  concrete 
dock.  This  Is  referred  to  as  a  smeared  reinforcement  model.  No 
consideration  was  given  to  modeling  various  discontinuities  (l.e., 
stairwells,  etc.)  or  the  shear/construction  joints  of  the  drydock. 

The  drydock  was  broken  Into  seven  substructures.  These  were  the 
right  main  sidewall,  the  left  main  sidewall,  the,  right  end  closure,  the 
left  end  closure,  the  pumphouse,  the  head  end,  and  the  calssori.  The 
level  of  discretization  of  the  substructure  was  decided  upon  using  the 
criteria  that  the  aspect  ratio  of  any  element  could  not  exceed  4.0  to 


pro«ot<*  /»f!pqunt<»  ncciirflr.y  In  rnmputoH  nn<l  the  minimum  number 

of  solid  elements  thtoiinh  fl  cross  sei  f  lo  t  uns  to  be  four  so  that  bending 
moments  could  he  mlnlmully  enptured.  The  computet  lone  1  s(7.n  of  the  sub¬ 
structures  wes  bfllanced  rtgrtlhst  the  computet  lone  1  size  of  the  main  system 
model  by  keeping  the  dof  nnd  brfndwidth  of  the  n  In  model  comparflble.  to 
the  dof  end  bnndwtdth  of  nny  substructure.  F.xei,  ties  of  the  resulting 
substructure  me ’els  are  presented  in  Figures  Wfn)  through  17(f). 

All  of  the  substructures  were  crented  on  n  KompnterVlslon ' (Prime 
Computer,  Inc.,  1986)  CAD/CAM  system  using  finite  element  modeling  soft¬ 
ware.  This  software  was  chosen  because  of  its  familiarity,  ease  of 
translation  of  data  files  for  use  witli  ADINA  (the  analysis  program  to  be 
used),  and  it  had  graphics  suitable  for  presentation  requirements.  How¬ 
ever,  this  CAD/CAM  system  proved  to  be  inadequate  for  large-scale  modeling. 
It  had  two' serious  shortcomings.  First,  execution  of  the  software  was 
very  slow.  Predictably,  as  the  number  of  sol  id  elements  in  the  drydock 
model  Increased  the  speed  decreased.  However,  the  system  required  2  hours 
to  generate  four  repeated  slices  of  250  elements  in  a  sidewall  substructure, 
and  it  required  up  to  10. minutes  to  renumber  200  nodes.  Second,  the 
capability  of  making  parametric  changes  during  construction  of  the  sub¬ 
structure  models  was  not  supported  by  this  system  Data  for  the  nodes 
and  foundation  spring  elements  (appr.oximately  1,000  each)  had  to  be  input 
one  at  a  time.  The  times  required  to  finish  some  of  the  substructures 
are  listed  in  Table  3. 


Table  3.  Time  Invested  in  Drydock  Substructure 
Model  Construction 


Substructure 

.Time  to  Complete  ■ 
(weeks ) 

Right  and  left  main  sidewall 

5 

Right  and  left  end  closure 

12 

Caisson 

6 

Head  end  (not  completed) 

4 

Pump  house  (not  started) 
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Figure  17.  Continued. 


RgurelZ.  Continued. 
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Rgure  17.  Continued. 


It  was  later  decided  to  switch  the  analysis  package  used  from  ADINA 
to  GIFTS.  This  was  done  because  ADINA  did  not  have  a  usable  substructur¬ 


ing  capablHty  that  suited  our  problem.  To  use  GIFTS,  the  model  data 
files  had  to  be  translated  from  ComptiterVlslon  to  GIFTS  format.  A  FORTRAN 
program  was  written  to  do  this.  It  required  8  weeks  to  write  the  file 
translation  computer  code,  translate,  the  flies  into  GIFTS  format,  and 
verify  the  models.  This  was  all  done  on  the  PC  version  of  GIFTS.  Another 
6  weeks  was  devoted  to  using  the  GIFTS  equation  solver  on  each  set  of 
substructure  equations  to  verify  that  the  substructure  model  was  valid. 

An  attempt  was  made  to  eliminate  the  Internal  nodes  of  the  right 
main  sidewall  substructure  using  GIFTS.  However,  after  a  significant 
effort.  It  was  found  that  the  GIFTS  substructure  reduction  capability 
could  not  cope  with  constrained  nodes  or  springs  on  the  boundaries  of 
the  sidewall.  Analyzing  the  drydock  using  substructure  modeling  with 
GIFTS  became  Impossible  because  it  could  not  model  the  constrained  nodes 
In  the  main  model  without  assigning  the  main  model  all  the  dof  of  the 
substructures  combined.  This  would  defeat  the  idea  of  substructure 
modeling.  Then  It  was  decided  to  at  least  make  a  simpler  "free-free" 
(l.e.,  no  restraining  nodes  and  no  soils  modeled)  model  for  purposes  of 
mode  shape  analysis  using  the  substructure  technique.  However,  problems 
were  found  In  the  GIFTS  Implementation  of  the  substructure  reduction 
algorithm  for  eigenvalue  analysis.  They  were  documented  and  forwarded 
to  GIFTS  support  engineers. 

A  frequency  analysis  of  the  sidewall  modeled  as  an  independent 
structure,  with  soil  foundation  and  boundary  constraints  included,  was 
performed  for  a  single  one-element  slice  of  the  right  main  sidewall  sub¬ 
structure  model.  Two  versions  of  the  model  were  employed  (see  Figure 
9);  one  used  solid  elements  to  represent  the  soil  and  the  other  used 
spring  elements  to  represent  the  soil.  The  eigenvalue  analysis  was  com¬ 
pleted  on  these  two  models  using  GIFTS.  (The  resulting  natural  fre¬ 
quencies  and  mode  shapes  were  presented  in  the  previous  section  of  this 
report  where  two-dimensional  models  were  discussed,  see  Figure  10.) 

Large-scale  analysis  of  a  fully  three-dimensional  solid  element 
model  of  the  drydock  using  the  substructure  concept  was  not  possible, 
due  to  the  various  problems  and  limitations  with  commercial  software 
products  mentioned  above  Involving  substructure  modeling,  instead,  a 
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much  simplified  three-dimensional  model  of  the  drydock  as  an  independent 
structure  was  constructed  using  plate  bending  finite  elements  exclusively. 
Then  using  GIFTS,  the  eigenvalue  analysis  for  the  natural  frequencies 
and  natural  mode  shapes  for  this  three-dimensional  drydock  model  was 
easily  completed.  (The  results  are  presented  and  discussed  in  a  later 
section,  Drydock  Model  Study.) 

In  summary,  the  ComputerVlslon  CAD/CAM  system  did  not  provide 
sufficient  hardware  and  software  capability  for  construction  of  large, 
complex  three-dimensional  finite  element  models  within  any  reasonable 
length  of  time  and  cost.  The  specific  problems  are; 

•  Hardware  speed  is  too  slow. 

I  / 

•  Wireframe  and  solids  modeling  algorithms  are  Inadequate. 

•  Combinations  of  solid  and  structural  elements  cannot  be 
processed. 

•  Substructure  combinations  cannot  be  processed. 

It  is  anticipated  that  new  developments  in  the  solids  modeling 
Ind  istry  as  reflected  in  the  latest  version  of  the  PATRAN  Plus  (1988) 
computer  program  or  the  MOVIESTAR  computer  program  recently  developed  at 
Brigham  Young  University  under  Navy  sponsorship,  for  example,  will  over¬ 
come  the  above  problems  and  provide  the  necessary  software  capability  to. 
develop  fully  three-dimensional,,  substructured  models  of  Navy  drydocks. 

Had  the  three-dimensional  drydock  system  model  envisioned  here  been 
successfully  developed.,  it  would  have  possessed  several  hundred  thousand 
degrees  of  freedom.  It  is  apparent  that  commercial  finite  element  analysis 
software  will  haye  to  improve  substantially  if  the  substructure  theory 
approach  is  to  be  employed  to  manage  the  associated  large-scale  computa¬ 
tional  burden. 
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Caisson  Model  Study 

Various  finite  element  modeling  approaches  for  analyzing  the  steel 
caisson  of  drydock  No.  6,  Puget  Sound,  are  considered  and  demonstrated. 
Emphasis  is  given  to  caisson  seal  reaction  loads  assuming  mainly  that 
the  seal  Is  rigid. 

Equivalent  Flat  Plate  Finite  Element  Study  of  Caisson.  To  develop 
a  preliminary  understanding  of  the  behavior  of  the  caisson,  calculations 
were  made  on  an  equivalent  plate  which  is  simply  supported  on  three  sides. 
This  equivalent  plate  can  be  modeled  via  analytical!  methods  in  Timoshenko 
and  Woinosky-Krleger  (1959)  or  via  ntimerical  finite  element  methods.  In 
both  cases,  the  simply  supported  sides  and  the' hydrostatic  loads  can 
easily  be  defined.  The  analytical  model,  however,  Is  generally  limited 
to  rectangular  plates  and  assumes  classical  boundary  conditions  so  that 
the  caisson  seal  will  Incorrectly  carry  tension.  Alternatively,  the 
finite  element  model  allows  for  other  than  rectangular  shapes  (most 
caissons  are  trapezoidal)  and  more  sophisticated  boundary  conditions 
including  unilateral  contact  conditions  or  gapping.  Either  method, 
however,  requires  establishment  of  an  equivalent  plate  thickness.  The 
calculations  for  the  equivalent  plate  thickness  of  the  caisson  are  given 
In  Appendix  D.  The  result  used  for  the  study  was  determined  to  be  58.1 
inches . 

The  symmetric  half  of  the  (176-foot  4.5-lnch  by  62-foot  10.5-inch, 
width  versus  depth)  caisson  was  modeled  as  a  rectangular  plate.  Numerical 
calculations  v;ere  made  using  an  equivalent  flat  plate  finite  element 
model  of  the  caisson.  (The  analytical  flat  plate  solutions  were  also 
considered,  although  they  are  omitted  here  for  brevity.)  Initial  numeri¬ 
cal  tests  were  made  to  determine  convergence  rates  of  various  finite 
element  models  of  a  rectangular  plate  simply  supported  on  three  sides 
and  subjected  to  a  hydrostatic  load.  The  criteria  used  to  evaluate  the 
accuracy  of  the  models  is  the  centerline  deflection  at  the  top  of  the 
caisson  and  the  boundary  reaction  loads  on  the  caisson  seal.  The  boundary 
supports  are  assumed  to  be  rigid.  The.  four-node  shear  flexible  plate/shell 
elements  In  the  ABAQUS  finite  element  program  were  used. 


The  first  series  of  tests,  used  to' check  the  number  of  elements 
which  are  necessary  for  accurate  solution  convergence,  employ  the  fol¬ 
lowing  caisson  models:  (1)  model  1  uses  a  4  by  10  element  mesh,  (2) 
model  2  uses  a  28  by '20  mesh,  and  (3)  model  3  uses  a  42  by  30  mesh.  Key 
results  are  shown  in  Table  4  and  the  seal  responses  are  shown  in  Figures 
18  through  20,  respectively. 


Table  4.  Results  of  Finite  Element  Plate  Model' 
Convergence  Study 


Caisson 

Model 

No. 

Mesh 

Centerline 
Deflection 
(in. ) 

Seal  Loads, 

Top 
Corner 
(lb/ in.) 

Bottom, 
Corner 
(lb/ in.) 

Centerline 
' Bottom 
(Ib/ln.) 

1 

10  X  4 

1.633 

-4 1 . 066 

20,593 

-9,678 

2 

28  X  20 

1.633 

-62,277 

30,837 

-9,736 

3 

42  X  30 

1.633 

-77,203 

36,486 

-9,750 

It  is  noted  that  the  seal  response  shows  a  large  concentrated 
compression  at  thp  top  corner  and  a  large  concentrated  tension  at  the 
bottom  corner  (see  Figures  18  through  20).  The  seal  response  along  the 
vertical  side  shows  an  Increase  with  depth  as  expected  with  a  hydrostatic 
loading.  The  responses  along  the  bottom  show  a  very  slight  increase 
toward  the  centerline,  but  is  relatively  constant.  As  the  mesh  gets 
finer  the  responses  at  the  corners  tend  to  get  more  concentrated.  The 
finite  element  solution  is  converging  to  what  the  linear  theory  of  plates 
suggests,  an  infinite  concentrated  load  at  the  corner  of  the  free  and 
supported  edges.  The  centerline  deflection  and  response  of  the  seal 
along  the  bottom,  however,  show  little  change  with  mesh  size.  The  28  by 
20  mesh  was  chosen  for  further  analysis  based  upon  these  results. 

Caisson  model  4  is  a  28  by  20  symmetric  half  of  a  mesh  that  is 
trapezoidal  vice  rectangular.  The  width  across  the  top  and  the  depth 
are  unchanged,  but  the  width  across  the  bottom  is  smaller.  The  total 
hydrostatic  load  is  correspondingly  reduced  by  8  percent.  Results  are 
presented  in  Table  5  and  compared  with  those  from  model  2.  The  seal 
responses  are  shown  in  Figure  21. 
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Rgure  18.  Rectangular  Plate  10x4,  Seal  Reactions. 


centerline  deflection  =  1.633' 


Rgure  19.  Rectangular  Rate  28  x  20,  Seal  Reactions. 


77^  #/ln. 


A  nlight  lncre«»«  In  c«‘nt*•rlln^  cl»*f Icr.t Ion  In  noted  for  the  trape- 
zoldal  ahape."  The  compreanlvn  seal  loads  decrease  significantly  at  the 
top  corner  and  Increase  slightly  at  the  bottom  near  the  centerline. 
Also,  the  tension  at  the  bottom  corner  Increases  significantly.  Figure 
21  Indicates  that  although  the  magnitude  of  ttie  peaks  change,  the  basic 
behavior  is  the  same. 


Table  5.  Results  of  Plate  Model  Heometry  Study 


Seal  I>oads  > 

Caisson 

Model 

No. 

Mesh 

Center  1 ine 
Deflect  ion 
(in.) 

Top 

Corner 

(Ib/in.) 

Bottom 

Corner 

(Ib/in.) 

Center  1 i ne 
Bottom 
(Ib/in.) 

2 

28  X  20 
Rectangular 

1.6TT 

-62.277 

30,837 

-9,736 

4 

28  X  20 
Trapexo Ida  1 

1.701 

-40. ISO 

■>0,019 

-9,910 

Caisson  model  5  is  nonlinear  because  it  allows  gapping  (or  tension 
not  allowed)  along  the  seal.  The  results  are  shown  in  Figure  22.  They 
Illustrate  how  gaps  may  tend  to  form  between  the  caisson  and  seal  near 
the  bottom  corner  (this  area  was  in  tension  In  previous  models).  Further, 
the  compressive  load  acting  on  the  seal  is  more  concentrated.  The  length 
of  seal  that  has  released  from  the  support  is  Ik  feet  of  the  total  142-foot 
width  with  a  maximum  gap  of  0.276  inches  at  the  corner.  Additional  results 
are  given  in  Table  6. 

In  s'lmmary,  the  nonlinear  caisson  model  shows  an  increase  in  maximum 
deflection  and  an  Increase  in  compressive  loads  on  the  seal.  The  maxi¬ 
mum  compression  occur Ing  along  the  bottom  is  of  particular  Interest.  In 
the  linear  model,  the  maximum  compression  along  the  bottom  occurs  at  the 
centerline.  The  nonlinear  model  results  show  the  maximtim  compression 
occurring  near  the  separation  point.  Similar  behavior  is  seen  along  the 
vertical  edge  of  the  seal. 
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Table  6.  Results  of  Linear  Model  Versus 
Nonlinear  Gapping  Model 


Seal 

Loads 

Caisson 

Model 

No. 

Analysis 

Centerline 
Deflection • 
(In.) 

Top 
Corner 
(lb/ In.) 

,  Maximum 
Bottom 
(Ib/ln.) 

4 

Trapezoidal 

linear 

1.701 

-40,159 

-9,910 

5 

Trapezoidal 

nonlinear 

2. 143 

-49,603 

-16,397. 

Maximum  for  the  nonlinear  run  does  not  qccur  at 
the  centerline  (see  Figure  22). 


Three-Dimensional  Finite  F.lement  Model  Study  of  Caisson.  The  sym¬ 
metric  half  of  the  caisson  for  drydock  No.  6,  Puget  Sound,  was  modeled 
using  detailed  three-dimensional  finite  element  technology.  The  internal 
frame  was  modeled  beam  for  beam  as  shown  In  Figure  23.  The  plate  struc¬ 
ture  consists  of  the  external  plates  that  are  stiffened  with  small  angle 
members  and  Internal  plates  (bulkheads)  that  are  used  for  stiffening  the 
caisson  and  separating  compartments,  and  these  are  modeled  as  shown  in 
Figure  24.  The  external  plates  are  modeled  approximately  with  plates  of 
equivalent  thickness.  The  caisson  Is  supported  along  the  outside  edge 
of  the  seal, and  Is  subjected  to  gravity  and  hydrostatic  loads. 

Model  6  Incorporates  the  Internal  frame  and  shell  structure  of  the 
caisson  as  described  and  also  includes  the  stiffness  of  the  concrete 
placed  In  the  bottom  of  the  caisson  for  ballast.  The  computed  seal  loads 
are  shown  In  Figure  25.  The  concrete  stiffness  removes  the  compression 
spike  near  the  bottom  corner,  but  increases  the  tension  carried  at  the 
corner.  Key  results  of  model  6  are  compared  with  the  equivalent  plate 
model  trapezoidal  (model  4)  in  Table  7 

Model  6  shows  a  significant  reduction  in  deflection  and  in  the 
seal  loads.  The  maximum  compression  along  the  bottom  Is  nearly 
unchanged,  but  the  location  predicte<l  by  the  equivalent  plate  model  was 
at  the  centerline  whereas  model  6  predicted  a  location  between  the 
bottom  corner  and  the  centerline. 
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Figure  24.  Plate  Structure  of  Caisson  (Symmetric  Hatf). 


16,080  #/in. 


Rgure  25.  Seal  Reactions  •  Linear  Calculation. 
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Table  7.  Results  Comparing  Equivalent  Plate  and  3-D 
Finite  Element  Models  of  Caisson 


Caisson 

Model 

No. 

‘  Mode 1 

Type 

Centerline 
Deflection 
(In.  ) 

Seal  Loads 

Top 
Corner 
(Ib/ln. ) 

Bottom 

Corner 

(Ib/ln.) 

Maximum 

Bottom 

(Ib/ln.) 

4 

Equivalent 

1.701 

-40,159 

50,019 

-9,910 

Plate 

6 

3-D 

0.643 

-25,956 

16,080 

-9,870 

Maximum  for  the  nonlinear  run  does  not  occur  at  the 
centerline  (see  Figure  25). 


Model  7  is  a  nonlinear  analysis  version  of  model  6  lhat  allows 
gapping.  The  seal  response  is  shown  in  Figure  26.  Results  from  the 
nonlinear  equivalent  plate  trapexoidal  model  (model  5)  are  compared  to 
model  7  In  Table  8. 


Table  8.  Results  Comparing  Nonlinear  Equivalent  Plate  and 
Nonlinear  3-0  Finite  Element  Models  of  Caisson 


■ 

Seal  Loads 

Caisson 

Model 

No. 

Ana lys is 

Centerline 
Deflection 
(  in.  ) 

Maximum 
Gap 
(in.  ) 

Top 

Corner 

(Ib/in.) 

Maximum 

Bottom 

(Ib/in.) 

5 

Equivalent 

Plate 

2.  143 

-0.276 

-49,603 

-16,397 

7 

3-D 

0.763 

-0. 102 

-31,381 

-10,713 

s 

Maximum  for  the  nonlinear  run  does  not  occur  at  the 
centerline  (see  Figure  26). 
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Model  8  estimates  the  behavior  of  the  caisson  bearing  against  a 
flexible  seal.  The  flexible  seal  Is  implemented  In  a  crude  fashion  via 
linear  springs.  The  results  from  model  8  show  that  a  flexible  seal  dis¬ 
tributes  the  load  uniformly.  Further,  there  Is  no  gapping.  The  center- 
line  deflection  Is  2.4A5  Inches.  Deflections  at  the  top  corner  and  the 
bottom  centerline  are  1.451  and  0.669  inches,  respectively.  These  results 
are  shown  In  Figure  27. 

Natural  Vibration  Analysis  of  a  Caisson 

A  natural  vibration  analysis  considering  the  symmetric  half  model 
of  the  caisson  (Figures  23  and  24)  was  performed.  The  natural  frequencies 
and  periods  for  the  first  five  symmetric  modes  are  given  in  Table  9. 

The  corresponding  natural  mode  shapes  are  shown  In  Figures  28(a)  through 
28(e). 

Table  9.  Caisson  Natural  Frequencies  and  Periods 


Mode 

No. 

Frequency 

(Hz) 

Period 

(sec) 

1 

9.041 

0.1106 

2 

20.579 

0.0486 

3 

24.904 

0.0402 

4 

24.911 

0.0401 

5 

24.918 

0.0401 

The  first  mode  shape  of  vibration  (Figure  28a)  is  similar  to  the 
shape  the  caisson  would  take  when  subjected  to  a  hydrostatic  load.  It 
would  be  expected  to  bn  the  dominant  mode  in  the  dynamic  response  to  an 
earthquake  load  or  to  a  uniform  blast,  load.  Tlie  second  and  third  modes 
(Figures, 28b  and  28c)  are  similar  In  shape  in  the  vertical  direction  to 
higher  modes  of  a  vertical  cantilever  wall.  Tlie  fourth  mode  (Figure  28d) 
is  a  local  deformation  mode  In  the  hull,  as  is  the  fifth  mode  (Figure 
28e)  which  also  contains  a  global  horizontal  deformation  mode. 


Rgure  27.  Seal  Reactions,  Nonlinear  Calculation  -  Rexible  Seal. 
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Drydock  Model  Study 


A  simplified  three-dimens ionnl  finite  element  model  of  drydock  No.  6 
at  Puget  Sound  was  constructed  to  demonstrate  the  fundamental  nature  of 
the  natural  modes  of  deformation  inherent  in  the  dynamic  response  of  a 
drydock.  One-  and  two-dimensional  models,  as  discussed  in  previous 
sections  of  this  report,  yielded  limited  information  in  this  respect. 
Despite  the  inability  to  effectively  construct  a  fully  three-dimensional 
model  (with  substructure-level  detail)  in  a  reasonable  length  of  time,, a 
simplified  model  was  constructed  which  was  sufficient  to  demonstrate  the 
natural  f»'equencies  an(.  ';.;t  iral  mode  shapes  of  a  large,  thin-walled 
drydock. 

Plate  Bending  Finite  Element  Model  of  a  Drydock.  The  three- 
dimensional  model  of  the  drydock  is  hased'on  the  cross-sectional  ideal¬ 
ization  of  drydock  6  (see  Figure  1)  used  to  construct  the  one-dimensional 
continuous  beam  model  of  ttie  drydock,  This  latter  model  was  discussed 
in  the  first  section  of  this  report  where  its  natural  frequencies  and 
mode  shapes  were  presented;  these  data  will,  therefore,  be  comparable 
with  the  calculated  values  of  natural  frequencies  and  mode  shapes  for 
the  present  three-dimens ion  il  model. 

The  present  model  was  developed  and  analyzed  using  the  GIFT.S 
general  purpose  finite  element  program.  The  model  is  shown  in  Figure  29 
and  consists  of  over  780  four-node  quadrilateral  plate  bending  finite 
elements.  A  summary  of  the  model  is  presented  in  Table  10.  These  data 
have  been  described  and  used  in  defining  models  previously  discussed  for 
the  concrete  drydock  and  for  the  steel  caisson.  '  The  three-dimensional 
model  is  consistent  with  the.  one-dimensional  continuous  beam  model  with 
tip  masses  that  represent  the  inertia  of  the  endwa 1 1  and  caisson.  The 
caisson  section  of  the  mode]  is  consistent  with  the  plate,  bending  finite 
element  model  used  to  study  door  seal  loads.  The  overall  stiffness 
distri’ ut ion  and  mass  distribution  of  drydock  No.  6  is  reasonably  well 
replicated  with  this  simple  three-dimensional  strncttiral  idealization. 


Rgure  29.  Three-Dimensional  Platfe  Bending  Finite  Bement  Model 
of  Drydock  6  at  Puget  Sound  NSY. 


A  .similar,  simple  Idealization  of  the  surrounding  soil  and  founda¬ 
tion  of  the  drydock  Is  omitted  from  the  three-dimensional  model  for  the 
sake  of  expediency.  From  the  previous  study  of  one-dimensional  models, 
it  was  shown  that  the  natirfll  node  shapes  (not  natural  frequencies)  of 
the  model  were  Independent  of  the  soil  foundation.  It  may  be  expected 
that  an  approximately  similar  situation  exists  for  the  present  three- 
dimensional  model. 


Natural  Frequuncies  and  Natural  No<l«*  Shapes.  The  first  10  nonzero 
natural , frequencies  and  natural  mode  shapes  are  shown  in  Figures  30(a) 
through  30(j).  Each  figure,  includes  four  views  of  the  mode  shape  infor¬ 
mation  along  with  the  value  of  the  natural  frequency.  These  different 
views  aid  in  determining  the  nature  of  the  deformation  in  each  mode. 
Recall  that  bending  moment  magnitude  is  proportional  to  magnitude  of 
curvature  in  the  walls  and  floor  of  the  deformerl  drydock.  Thus,  the 
location  and  relative  intensity  and  direction  of  bending  moment  in  the 
drydock  walls  and  floor  associated  with  each  mode  shape  can  be  detected 
by  scanning  the  deformed  finite  element  mesh  for  areas  of  high  curvature. 

The  main  observation  is  that  these  natural  mode  shapes  are  primarily 
composed  of  longitudinal  deformation  modes.  These  modes  are  somewhat 
similar  to  the  modes  obtained  from  the  one-dimensional  continuous  beam 
model,  but  they  are  entirely  different  from  those  obtained  from  the.  two- 
dimensional,  plane  strain  finite  element  model. 
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Figure  30.  Condnued. 
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Figure  30.  Continued. 


The  first,  and  fundnmentfl  1 ,  of  vibration  fov  the  drydock  model  , 

occurs  at  a  frequency  of  0.3fl  hertz,  which  is  very  low.  Notice  that  all 
the  generators  in  the  deformed  mesh  are,  straiglit  lines  Indicating  that 
very  little  or  no  bending  occurs  in  the  drydock  for  this  mode.  This 
mode  is  a  longitudinal  torsional  mode  in  whiclj  sliear  forces  predominate. 
Recall  that  the  one-dimensional  continuous  beam  model  Is  Incapable  of 
predicting  longitudinal  torsional  natural  modes  of  vibration. 

The  second  mode  shape  is  a  longitudinal  bending  mode  which  is  very 
similar  to  the  first  fiindamental  mode  shape  calculated  using  the  one- 
dimensional  'ontlnuoiis  beam, model.  •  It  is  of  interest,  therefore,  to 
also  compare  their  corresponding  natural  frequencies,  0.87  hertz  versus 
1104  hertz.  The  former  value  from  the  finite  element  model  should  be 
more  accurate.'  Thus,  the  one-dimensional  model  prediction  of  this 
natural  frequency  is  about  20  percent  high. 

The  higher  modes  are  increasingly  more  localized  combinations  of 
longitudinal  and  crOss-sect iona 1  bending'  modes . 

TnREF.-DTMF.NSIONAr.  NONT.TNEAR  FINTTR  F.IJilMENT  ANALYSTS 
OF  A  PIF.R  DECK  SCAUR  MODEL 

Finite  element  technology  can  provide  structural  assessment  of  other 
waterfront  facilities  that  impact  fleet  reafliness,  such  as  piers  and 
wharfs.  The  computer  code  ADIN'.  (ADTNA  Rngineerieg,  Tnc;,  1987)  is  used 
to  predict  the  response,  including  load-deflection,  stress  distribution, 
extent  of  cracking,  and  failure,  mode,  of  a  1/3-scale  reinforced  concrete 
pier  deck  model.  A  three-dimensional  nonlinear  reinforced  concrete  model 
is  reqtiired  to  accurately  represent  the  progressive  crack  formation  during 
loading,  and  to  reflect  the  consequent  stiffness  degradation. 

'  Agreement  with  experimental  data  is  essential  to  the  assessment  of 
accuracy  of  analytical  predictions..  Fortunately,  extensive  data  are 
available  for  a  i/3-scale  model  of  a  Navy  reinforced  concrete  pier  deck 
recently  built  and  tested  to  failure  in  related  exploratory  research  in 
structural  assessment  at  NCF.I..  In  the  following,  these  data  ate  used  to 
assess  the  accuiacy  of  three-dimensional  nonlinear  finite  element  tech¬ 
nology  for  reinforced  concrete  structures. 
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Pier  D«ck  Scale  Model 

Description.  The  model  has  five  spans  supported  by  six  pile  bents, 
and  is  showi.  in  Figure  31.  Dimensions  and  reinforcement  used  are  detailed 
in  Figure  32.  The  deck  depth  is  5.375  Inches,  representing  a  16-lnch 
prototype  depth.  The  main  reinforcement  consists  of  #3  grade  bars  and 
D5  deformed  wire,  representing  #9  and  //6  bars  In  the  prototype  pier. 

The  pier  is  designed  for  one-way  action,  hence  the  main  reinforcing  Is 
longitudinal.  Secondary  transverse  steel  is  also  included  at  top  and 
bottom. 

Material  Properties.  The  concrete  mix  is  indicated  In  Table  11. 

The  maximum  aggregate  size  was  scaled  down  to  3/8  inch.  A  higher 

relative  content  of  lime  aggregate  and  the  use  of  a  high-range  water- 

reducing  admixture  (superplasticize.r)  allowed  the.  concrete  to  be  pumped. 

Concrete  and  steel  properties  are  indicated  in  Table  12.  The  maximum 

concrete  compressive  stress  of  7,700  psi  was  reached  at  a  strain  of 

0.0027.  The  content  of  fine  measured  modulus  of  elasticity  was  found  to 

be  slightly  lower  than  the  traditional  value  of  57,000  /T  for  normal 

c 

weight  concretes.  The  measured  Poisson  ratio  was  0.15  (typical  values 
range  from  0.15  to  0.22,  ASCE,  1982). 

Test  Setup.  A  single  centered  patch  load  was  applied  to  the  center 
span.  Patch  dimensions  were  8  by  8  inches,  simulating  a  24-  by  24-lnch 
crane  outrigger  load  in  the  prototype  pier.  Tlie  load  was  applied  in 
four  steps,  up  to  30,  60,  and  90  kips,  then  to  failure.  At  each  of  the 
first  three  steps,  the  specimen  was  loaded  then  unloaded  ten  times.  The 
load  was  applied  through  a  100-ton,  hollow  ram,,  hydraulic  jack  pushing 
on  a  8-  by  8-  by  1-1/2-inch  steel  plate,  and  reacting  against  a  high- 
strength  steel  rod  anchored  to  the  laboratory  floor.  The  rod  had  been 
previously  calibrated  to  150  kips. 

Deflections  were  measured  at  the  center  of  the  load  patch  using 
linear  variable  differential  transformers  (LVDTs). 
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Table  11.  Concrete  Mix* 


Material 

Weight  Per 
Cubic  Yard 

Type  I  Cement 

658  Ib/y^ 

Fly  Ash 

100  Ib/y^ 

Water 

350  Ib/y^ 

3/8"  Gravel 

1,080  Ib/y^ 

Sand 

1,740  Ib/y^ 

*Water  reducing  admixture 
(Sikament  86)  per  manufacturer. 


SR-4  paper-backed  concrete  strain  gages  were  epoxied  to  the  com¬ 
pression  face  of  the  deck  slab  near  the  load  point  (1  Inch  away  from  the 
bearing  plate).  Weldable  strain  gages  were  attached  to  the  tension  rein¬ 
forcement  opposite  the  compression  gages.  These  strain  gage  pairs  provided 
a  measure  of  the  strain  field  across  the  deck  depth. 


96 


Nuaerlcal  Model 


Due  to  the  existing  symmetry,  only  one  quarter  of  the  physical  pier 
deck  model  is  modeled  numerically.  The  finite  element  mesh  used  is  shown 
in  Figures  33  and  3A.  Concrete  in  the  center  span  was  discretized  using 
five  layers  of  low-order  3-D  isoparametric  elements,  with  eight  nodes,  3 
degrees  of  freedom  per  node,  and  a  2x2x2  integration  order.  The  mesh 
used  is  shown  in  Figure  35. 

The  material  model  used  is  a  nonlinear  elastic  concrete  model  with 
compression  crushing  and  tensile  cracking  cut-off  with  strain  softening 
behavior.  Reinforcing  bars  are  represented  using  simple  truss  elements. 
This  method  is  referred  to  as  a  discrete  reinforcement  model.  Top  and 
bottom  reinforcement  nieshes  are  also  shown  in  Figure  35.  A  nonlinear 
elastic-plastic  material  model  was  used  for  the  steel. 

Contiguous  spans  were  modeled  using  linear  four-node  shell  elements 
with  5  degrees  of  freedom  per  node,  and  iincracked  concrete  properties. 

In  this  way  the  confinement  of  the  adjacent  spans  on  the  center  span 
could  be  adequately  represented  without  excessively  increasing  the 
computational  effort.  Transition  shell  elements  were  used  around  the 
center  span  to  connect  the  shell  to  the  3-D  elements.  The  model  had  a 
total  of  4,083  degrees  of  freedom. 

Owing  to  symmetry,  only  1/4  of  the  concentrated  vertical'  load  needs 
to  be  considered  in  the  analysis.  This  load  is  distributed  on  the  top 
corner  element,  which  has  dimensions  of  4x4x1  inch. 

The  pile  bents  rested  on  plywood  sheets,  which  were  also  Included 
in  the  finite  element  model  (see  Figure  33). 

Concrete  Con.stltutive  Laws  in  the  Post  Cracking  Range 

Reinforced  concrete  structural  failures  primarily  evolve  from  tensile 
cracking.  Thus  it  Is  Important  to  check  the  behavior  and  performance  of 
any  nonlinear  concrete  material  model,  especially  with  regard  to  its 
tensile  cracking  behavior,  before  the  model  is  employed  in  the  analysis 
of  a  structure. 
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Figure  34.  Finite  Element  Mesh,  Top  View. 


ELASTIC  PLASTIC  STEEL  ELEMENTS 


Figure  35.  Concrete  and  Steel  Meshes,  Center  Span. 
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PrellaJuary  Material  Model  Testing.  As  pnrt  of  the  preliminary 
computer  program  evaluation,  a  simple  tensile  Ltst  on  a  single  four-node 
2-D  concrete  element  was  carried  out.  The  element  was  pulled  vertically  ' 
frcm  the  top  two  nodes,' then  unloaded,  then  pulled  until  no  more  tensile 
stresses  were  transfer* ed.  All  other  degrees  of  freedom  were  deleted. 
From  Figure  36a  it  may  be  observed  that  the  unloading  branch  is  Incorrect 
because  it  is  not  directed  at  the  origin.  Since  theADINA  source  code 
is  provided,  it  was  possible  to  correct  this  deficiency  as  Figure  36b 
demonstrates.  The  >:aipa  corrections  were  made  to  the  3-D  concrete 
element. 

The  tensile  test  was  repeated,  this  time  allowing  for  lateral 
contraction  (Poisson  deformation)  by  releasing  two  lateral  displacement 
degrees  of  freedom.  No  solution  was  obtained  when  the  usual  energy  con¬ 
vergence  criterion  in  ADINA  was  used.  This  problei;^  was  circumvented  by 
switching  to  a  force  and  moment  convergence  criterion  in  ADINA. 

Constitutive  Laws.  For  concrete  in  monotonic  uniaxial  tension 
beyond  the  cracking  point,  ADINA  provides  a  linear  release  of  the' 
tensile  stresses  (Figure  36).  This  part  of  the  stress-strain  curve  is 
referred, to  as  strain  softening.  Related  basic  research  in  NCEL's 
Structural  Modeling  project  has  shown  that  a  nonlinear  tensile  stress 
release  is  more  realistic.  The  normalised  strain  softening  relationship 
used  is  shown  in  Figure  37a.  In  this  figure,  limit  strain  Is  the  strain 
beyond  which  no  tensile  stresses  are  transferred.  This  relationship  was 
Implemented  in  ADINA  as  an  option  (Malvar  and  Warren,  1989;  Malvar  and 
Fourney,  1990).  A  fracture  energy  of  0.S7  Ib/in.  (100  N/m)  Was  also 
assumed. 

Accepted  behavior  after  concrete  cracking  provides  for  shear 
stresses  to  be  transferred  across  the  crack.  However,  they  decay 
quickly’ with  increasing  crack  opening.  To  account  for  this  in  ADINA, 
the  shear  modulus  is  decreased  linearly  with  increasing  strain  perpen¬ 
dicular  to  the  crack,  down  to  a  minimum  value  which  remains  constant. 
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2D  CUBE  IN  TENSION 


2D  CUBE  IN  TENSION 


Figure  36.  Tensile  Behavior  Correction  of  ADINA  Concrete  Model. 
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NORMALIZED  SHEAR  MODULUS  CONCRETE  STRESS/TENSILE  STRENGTH 


(a)  Tensile  strain  softening  normalized  stress  versus  strain  law 


Small  Slip  Model  — ^  Input  to  ADINA 

(b)  Post  cracking  shear  modulus 


Figure  37.  Tensile  Strain  Softening  and  Post  Cracknig  Shear  Modulus. 


Experiment A 1  observations  show  that  ^he  shear  stiffness  decays  expo¬ 
nentially.  For  small  slip  values  in  .larticular,  the  decrease  in  shear 
modulus  G  can  , be  derived  from  Bazant  and  Gambarova  (1987): 


‘’nt  *  S  6- 


where:  0  ^  =  shear  stress  transferred  across  crack 
nt 


8^  ,=  crack  slip 


6^  =  crack  opening 


*2  =  *3 


-3  = 


T  =  0.245  f* 
o  c 


(N/tnm^) 


T  =  t  a  /Ca  +8 
u  o  o  on 


a  =  O.Ol  D 

o  a 


=  maximum  aggregate  size 


For  small  crack  openings,  this  can  be  simplified  to: 


a  ,  =  2.45 
nt  8 


(N/mm^) 


Further,  if  it  is  assumed  that  t  =  5  /s  and  =  6^/s, 

n  n  t  ,  t  ' 


where  s  is  a  prescribed  crack  spacing,  then  the  degradation  of  the  shear 
niodulus  with  crack  opening  is  governed  by: 


G  = 


nt  2.45 
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This  is  nquivalent  to  an  exponential  shear  stiffness  decay  with 
Increasing  crack  width,  as  shown  in  normalized  form  in  Figure  37b.  The 
bilinear  approximation  used  as  input  to  ADINA  attempts  to  match  the  rapid 
shear  stiffness  decay  at  small  crack  openings.  The  ratio  of  cracked 
shear  stif fness  to'the  initial  uncracked,  value  is  often  termed  the'shear 
retention  factor. 

Experiaental  Results 

Measured  ’lOad-Def lection  Behavior.  For  the  first  load-unload  cycle 
at  each  load  step,  the  ascending  part  of  load-deflection  curve  is  shown 
in  Figure  38.  The  slope  of  the  curve  (stiffness)  shows  a  sharp  decrease 
at  around  30  kips,  representative  of  extensive  concrete  cracking.  Each 
successive  loading  or  unloading  is  linear  if  the  load  is  smaller  than 
the  maximum  load  previously  attained.  If  this  maximum  load  is  exceeded, 
the  load-deflection  path  shows  a  further  decrease  in  slope  and  follows 
an  envelope  similar  to  a  monotonlc  loading  test. 

In  the  final  cycle  to  failure,  a  maximum  load  of  121  kips  was 
obtained,  at  a  deflec^on  of  0.32  inch. 

Observed  Crack  Patterns  and  Failure  Node.  Crack  maps  were  recorded 
at  different  load, levels,  and  at  failure.  Cracks  at  the  top  and  bottom 
of  the  deck  are  reported  in  Figures  39  and  40,  respectively.  Top  cracks 
are  almost  circular  and  concentric  around  the  load  point.  Bottom  cracks 
are  radial  until  punching  failure  takes  place. 

Brittle  punching  shear  failure  occurred, at  121  kips.  On  the  deck 
top,  section  cracks  developed  just  around  the  load  footprint  at  an  Initial 
angle  of  about  45  degrees  with  the  vertical.  Figure  40  shows  the  inter¬ 
section  of  these  cracks  with  the  deck  bottom.  Deflections  in  the  slab 
at  failure  were  not  high  enough  to  develop  membrane  forces. 
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Figure  38.  Lx>ad  Deflection  History  at  Center. 
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SUPPORT  EDGE 


SPAN  3  -  TOP 


■  CYCLE  15,  30  KIPS 
—  CYCLE  I  90  KIPS 
—  CYCLE  20,  90  KIPS 
—  FAILURE 

Figure  39.  Experimental  Crack  Patterns  on  Deck  Top. 
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SPAN  3-  BOTTOM 
*  CYCLE  15.  30  KIPS 

—  CYCLE  1.  90  KIPS 

—  CYCLE  20,  90  KIPS 

—  FAILURE 


Figure  40.  Experimental  Crack  Patterns  on  Deck  Bottom. 


Measured  Concrete  and  Steel  Strain.  For  the  first  cycles  to  30, 

60,  and  90  kips,  as  well  as  for  the  cycle  to  failure,  the  strain  in  the 
extreme  concrete  compression  fiber  (deck. top)  and  the  strain  in  the 
tension  steel  next  to  the  load  point  were  recorded  (Figures  41  and  42, 
respectively).  Values  for  both  strains  kept  close  to  each  other  during 
testing,  indicating  that  the  neutral  axis  in  tliat  section  remained 
approximately  equidistant  from  the  deck  top  and  the  tension  steel. 

Numerical  Results 

Predicted  I«oad-neflection  Behavior.  In  the  predicted  load-deflection 
graph  (Figure  35),  a  marked  variation  in  the  slope  at  load  step  2  reflects 
a  prediction  of  extensive  cracking  of  the  concrete  elements.  After  load 
step  6,  the  numerical  response  is  totally  dependent  on  the  shear  retention 
factor  (i.e.,  the  prescribed  degradation  of  the  shear  modulus,  G) .  When 
the  shear  retention  factor  was  set  to  zero,  brittle  failure  was  predicted 
at  load  step  6.  When  post  cracking  shear  retention  was  introduced  (accord 
ing  to  Figure  37),  the  remainder  of  the,  curve  was  obtained  with  a  brittle 
failure  mode  predicted  at  a  load  of  122  kips  (load  step  9),  which  is 
within  1  percent  of  the  measured  failure  load. 

Predicted  Crack  Patterns  and  Failure  Mode.  The  numerical  crack 
maps  at  vhe  deck  top  and  bottom  are  shown  in  Figures  43  through  45. 

Crack  surfaces  and  their  orientation  at  integration  points  are  repre¬ 
sented  by  two  parallel  circles.  On  the  deck  top  (Figures  43  and  44), 
horizontal  in-plane  cracks  indicate  delamination  of  the  concrete  com¬ 
pression  zone  near  the  load  due  to  shear  stresses.  Away  from  the  load, 
vertical  cracks  due  to  flexural  stresses  form  almost  circular  patterns 
about  the  load.  On  the  bottom,  radial  surface  cracks  eminate  from  the 
load  point  outward  prior  to  failure,  as  'bown  in  Figure  45. 

At  load  step  6,  Figure  46  Indicates  the  initiation  of  a  conical 
section  crack  around  the  load,  typical  of  the  punching  shear  failure 
observed  experimentally  (the  load  is  applied  on  the  top  right  element  in 
Figure  46). 
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LOAD  (KIPS) 


First  Cycle  30  kips  — First  Cycle  60  kips  First  Cycle  90  kips 
Cycle  to  Failure  FEM  Concrete  Strain 


Figure  41 .  Measured  Versus  Computed  Concrete  Strain  History. 
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Figure  43.  Predicted  Crack  Patterns  on  Deck  Top  at  90  Kips. 
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Figure  44.  Predicted  Crack  Patterns  on  Deck  Top  at  Failure. 
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Figure  45.  Predicted  Crack  Patterns  on  Deck  Bottom  at  Failure. 
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Figure  46.  Punching  Shear  Cracks. 


Predicted  Concrete  and  Steel  Strain.  The  computed  concrete  strain 
on  the  deck  top,  near  the  load  point,  and  for  monotonic  loading.  Is  com¬ 
pared  with  measured  strain  (Figure  41,  FF.M  concrete  strain).  At  step  6, 
quasl-horlzontal  cracking  occurred  in  this  location  (upper  right  element 
In  Figure  43),  with  corresponding  large  element  deformations.  Figure  42 
shows  the  tension  steel  strain  for  the  same  case.  Prior  to  concrete 
cracking  In  tension,  the  steel  deformation  is  constrained. 

The  three-dimensional  finite  element  model  also  allows  for  convenient 
evaluation  of  stress  and  strain  distributions  through  the  thickness  of 
the  pier  deck,  as  Indicated  in  Figure  47  for  the  pier’  deck  center.  The 
concrete  compression  stress  is  shown  to  decrease  almost  linearly  down  to 
the  neutral  axis,  below  which  tension  stresses  and  cracking  are  apparent. 
The  steel  stress  Is  also  Indicated. 


REINFORCED  CONCRETE  PIER  DECK 


Figure  47.  Computed  Stress  Profile  at  Center. 
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Discussion 


Iioad-Dcflectlon  Behavior.  Tlie  computed  precracking  branch  (up  to 
30  kips)  follows  the  experimental  results  very  closely.  Beyond  cracking, 
the  finite  element  model  is  consistently  stiffer  (see  Figure  38).  Three 
causes  may  account  for  this  discrepancy.  First,  the  limited  amount  of 
stiffness  reformations  performed  around  the  cracking  point  may  not  be 
sufficient  to  represent  all  the  cracking  that  occurs  simultaneously. 

The  prescribed  relative  closeness  pf  load  steps  1,  2,  and  3  is  an  attempt 
at  preventing  this'  problem.  Second,  the  low-order  isoparametric  elements 
used  typically  yield  a  somewhat  stiffer  response.  Third,  if  the  assumed 
fracture  energy  is  larger  than  the  actual  value,  the  response  beyond  the 
cracking  point  will  be  stiffer. 

At  a  load  of  90  kips  (load  step  6),  the  finite  element  model  shows 
a  strong  softening  of  the  response.  Tlie  experimental  results  show  a 
rapid  stiffness  degradation  under  cycling  at  that  same  load. 

At  ultimate,  the  numerical  solution  did  not  converge,  and  the  model 
was  unable  to  represent  the  deformation  at  con.stant  load  which  preceded 
brittle  failure.  The  failure  load  was,  however,  accurately  predicted  by 
the  finite  element  model. 

Crack  Patterns  and  Failure  Mode.  For  the  deck  top,  the  crack  for* 
mation  in  the  finite  element  model  simulates  very  closely  the  actual 
patterns  (compare  Figures  43  and  44  with  Figure  39).  The  higher  numeri¬ 
cal  crack  density  is  due  to,  the  assumption  of  perfect  bond  between  con¬ 
crete  and  steel  reinforcement.  Tf  bond-slip  effects  were  to  be  simulated, 
these  cracks  would  merge  into  fewer,  wider  cracks,  which  would  be  more 
realistic.  Coalescence  of  the  cracks  when  considering  bond-slip  effects 
has  been  demonstrated  (Rots,  1985;  RIT.F.M  Technical  Committee,  1987). 

The  finite  element  model  also  indicates  the  presence  of  horizontal 
cracks  under  the  load,  due  to  the  high  horizontal  compression  in  the 
concrete.  These  cracks  could  cause  failure  by  de lamination  of  the  con¬ 
crete  cover. 

The  finite  element  model  was  able  to  successfully,  predict  the  type 
of  failure,  and  to  show  the  effects  of  shear  retention  on  the  pier  deck 
response. 


Concrete  and  Steel  Strain.  Experimental  concrete  strain  readings 
on  the  deck  top  (Figure  41)  show  a  progressive  softening  with  successive 
load  cycles.  The  finite  element  results  shj)w  a  similar,  though  st if fer 
behavior  in  the  first  few  steps,  and  then  a  sudden  softening  after  step 
6.  The  excess  stiffness  may  be  attributed  to  the  use  of  low-order  iso¬ 
parametric  elements,  and  to  a  reduced  number  of  stiffness  reformations 
around  cracking.  The  sudden  softening  after  step  6  indicates  the  forma¬ 
tion  of  a  horizontal  crack  beneath  the  loafi,  representative  of  delamination 
in  the  compression  zone.  This  is  not  apparent  in  the  test  data.  It  may 
be  an  artifice  due  in  part  to  the  co.irse  discretization  around  the  load 
patch,  where  the  entire  load  is  distributed  over  the  surface  of  only  one 
element. 

The  finite  element  tension  steel  strain  prediction  forms  a  reasonably 
accurate  envelope  to  the  maximum  experimental  readings  (Figure  42).  The 
finite  element  morlet  sh<iws  a  constrained  deformation  prior  to  concrete 
cracking  in  tension  which  is  not  apparent  in  the  test  data.  This  may  be 
attributed  to  phenomena  which  were  hot  modeled  such  as  the  prior  existence 
of  tension  cracks  in  the  concrete  due  to  shrinkage,  and  to  some  low  load 
tests  (up  to  8  kips)  carried  out  initially. 
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SUMMARY  ANn  CONCIWSIONS 


A  sfit  of  Nflval  clvl 1 -striictural  onginoor ing  problems  Involving 
waterfront  facilities  was  considered  to  demonstrate  the  appropriateness 
of  advanced  three-dimensional  structural  analyses  for  the  statics  and 
dynamics  of  large  structiiral/geotechnical  systems.  A  drydock  is  a  pri¬ 
mary  example  of  such  a  facility.  There  is  a  continuing  Navy  requirement 
to  certify  ship  safety  while  it  is  in  drydock.  Thus,  a  drydock' system 
is  an  appropriate  subject  for  the  technoiogy  demonstration  contained  in 
this  report.  In  particular,  drydock  No.  6  at  Puget  Sound  NSY  was  selected 
because  it  is  a  relatively  modern,  fully  relieved  drydock  design  having 
thin  wall  and  thin  floor  sections.  Tt  is  therefore  a  more  suitable  sub¬ 
ject  for  advanced  analysis  methods  than  an  older  drydock  constructed 
mor.ol  i  thlca  1  ly. 

Natural  mode ■ shapes  and  natural  frequencies  are  fundamental  to 
structural  dynamics  of  engineering  systems.  Previous  analyses  of  Navy 
drydocks  were  not'  three-dimensional,  and  could  not,  therefore,  predict 
longitudinal  natural  modes  of  vibration.  Therefore,  a  fully  three- 
dimensional  finite  element  model  was  constructed  to  determine  the 
natural  modec-  ind  natural  frequencies  of  the  drydock.  The  results 
showed  tha*;  the^e  longitudinal  modes  are  most  important  because  they 
occur  -^it  the  liv?  end  of  the  frequenc5'  spectrum  for  drydocks,  and  are 
thus  more  apt  to  be  excited  by  typical  earthquake  ground  motion. 

Natural  mtifies,  calculated  using  n  two-dimensional  model  in  this  study, 
were  go'neraliy  higher  frequency  modes.  The  two-dimensional  model 
clearly  doss  rot  provide  the  same,  rirluiess  f>f  information  on  natural 
mode.s  of  dei'nr.iation  and  on  dynamic  response  of  drydocks.  Further,  the 
ciosevosp  of  u>wer  longitudinal  mode  frequencies  and  frequencies  typical 
of  ear'.hquake  spectra,  indicate  that  full,  three-dimensional  dynamic 
analysis  l.s  necessary  for  investigating  earthquake  response  of  drydocks. 

Since  establishing  ship  safety  is  a  paramount  goal  of  the  Navy's 
drydock  certification  program,  inclusion  of  a  simple  blocked  ship  model 
as  a  substructure  mode!  together  with  the  drydock  structural  model  was 
proposed  in  this  study.  A  .-Iniplif led  blocked  ship  model  was  developed 
and  used  to  predict  the  natural  mode  shapes  and  natural  frequencies  as 
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veil  as  the  dynamic  response  of  a  blocked  submarine  subjected  to  a 
simulated  floor  motion  excitation.  The  importance  of  knowing  the 
drydock  floor  motion  for  reliable  prediction  of  the  blocked  vessel's 
dynamic  response  was  emphasized. 

Measured  earthquake- Induced  drydock  floor  motion  data  available  to 
the  Navy  are  scant.  The  1987  Whlttler-Narrows  earthquake  excited  dry- 
dock  No.  1  at  Long  Beach  NSY  and  the  Indticerf  motion  was  measured. 

Analysis  of  the  data  Indicated  that  the  drydock  sustained  primarily 
rigid  body  motion,  and  therefore  It  would  not  have  been  exposed  to  any 
significant  structural  action.  Tlie  drydock  did  not  experience  a  strong- 
motion  earthquake  during  this  event.  Nonetheless,  the  blocking  system 
supporting  a  Naval  vessel  In  an  adjacent  drydock  did  experience  some 
damage.  In  general,  it  Is  expected  that  substantial  ground  motion  will 
Induce  both  dynamic  stresses  In  a  drydock,  and  motion  of  any  blocked 
vessels.  Earthquake  analysis  procedures  supporting  drydock  certifica¬ 
tion  should  therefore  consider  the  Interaction  between  the  dynamic 
response  of  blocked  vessels  and  the  dynamic  response  of  the  drydock. 

A  three-dimensional  nonlinear  finite  element  analysis  demonstration 
of  a  steel  drydock  caisson  substructure  has  provided  new  insight  Into 
the  behavior  of  forces  acting  on  the  door  seal  due  to  hydrostatic  loads. 
This  behavior  Is  different  from  that  which  tinderlles  current  caisson 
design  methods,  and  could  Improve  caisson  design  methods  and  maintenance 
procedures.  Further,  a  natural  vibration  analysis  using  the  three- 
dimensional  finite  element  model  of  the  caisson  was  conducted  to  suggest 
how  the  caisson's  vtilnernblllty  to  earthquake  and  blast  loads  could  be 
appropriately  Investigated  with  this  technology. 

A  three-dimensional  nonlinear  f in  i  te,  el ement  model  of  a  reinforced 
concrete  pier  deck  was  developed  and  employed  to  demonstrate  the  accuracy 
of  the  technology  in  predicting  the  response  of  the  structure  to  static 
load.  Concrete  and  steel  bars  were  modeled  using  a  nonlinear  concrete 
constitutive  model  and  a  discrete  reinforcement  model,  respectively. 

Model  predictions  for  the  following  detailed  responses  were  In  good 
agreement  with  experimental  measurements  and  observations:  monotonlc 
load-displacement  history;  progressive  crack  formation  and  final  crack 
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petternsj  progressive  stiffne-is  degradation  due  to  cracking;  a  punching 
shear  failure  mode;  and  concrete  and  steel  strain  histories  through  the 
deck  depth. 

Three-dimensional  linear  and  nonlinear  finite  element  technology  is 
appropriate  for  engineering  projects  alme<l  at  analysis  and  design,  and 
aimed  at  understanding  and  improving  serviceability,  vulnerability,  and 
survivability  of  drydccks  and  other  large  structural/geotechnical  systems 
along  the  waterfront. 
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A  BRIEF  REVIEW  OF  SOLIDS  MODELING  RELATED  TO 
3-D  FINITE  EIJC MENT  MODEL  CONSTRUCTION 

The  wireframe  geometric  model  Is  the  simplest  of  the  modeling  schemes 
and  requires  no  sophisticated  finite  element  preprocessing.  In  fact, 
determining  the  size  and  distribution  of  finite  elements  requires  much 
Ingenuity  on  the  part  Of  the  analyst.  An  example  of  a  wireframe  modeler 
Is  the  ComputerVlslon  system.  The  boundaries  of  the  geometric  model  are 
outlined  using  wireframes.  The  simplest  wireframe  element  consists  of 
two  nodes  connected  by  a  line,  scaled,  and  given  some  special  orientation. 
Two  or  more  nodes  can  occupy  the  same  point  In  .space  but  In  general  these 
wireframe  elements  have  no  connection  (or  connectivity)  with  other  wireframe 
elements.  The  amount  of  time  to  assemble  a  geometric  model  Is  greatest 
for  the  wireframe  method  since  each  line  has  to  be  Individually  Inserted. 
Shaded  Imaging  for  subsequent  visualization  requires  additional  work 
since  surfaces  also  have  to  be  defined  and  inserted. 

Development  of  finite  element  meshes  from  the  wireframe  model 
generally  requires  additional  work.  The  finite  element  modeler  generates 
meshes  using  primitive  shapes  (i.e.,  cubes,  boxes,  cylinders,  spheres, 
etc.).  The  wireframe  model  Is  subdivided  Into  these  primitive  shapes  to 
facilitate  generation  of  the  finite  element  mesh.  Once  the  shapes  or 
volumes  are  created,  the  mesh  Is  developed  by  defining  the  number  of 
elements  or  nodes  along  the  wireframe  edges  of  the  primitive  shapes. 

This  volume  Is  filled  with  elements  and  nodes  that  have  complete  con¬ 
nectivity  within  Itself  but  Is  otherwise  unrelated  to  other  volumes.  An 
adjacent  shape  Is  chosen  for  the  next  finite  element  discretization  but 
now.  In  addition.  It  requires  a  strategy  on  how  it  connects  to  the  previous 
volume.  The  number  of  nodes  (or  element  projections)  on  the  common 
surface  must  match  the  discretization  of  the  new  volume.  If  the  volumes 
are  different  sizes  (the  common  surface  is  only  part  of  the  surface  of 
the  larger  volume),  then  the  wireframe  model  must  be  further  modified 
(so  that  the  size  of  the  surfaces  matches)  or,  via  a  trick,  the  larger 
volume  must  be  discretized  so  that  the  nodes  match  up  at  the  common  sur¬ 
face.  The  meshing  Is  not  complete,  however,  because  the  two  volumes 
with  their  complete  connectivities  and  matching  common  surface  have  no 
connectivity  between  themselves.  A  further  step  of  merging  the  duplicate 
nodes  (nodes  occupying  the  same  point  in  space)  finally  establishes  a 
connectivity  between  the  elements  on  the  common  surface  and;  hence,  the 
volumes . 

A  second  scheme  Is  called  the  surface  geometric  model.  It  Is  a 
"solld-like"  model  because  It  uses  solid  primitives  for  assembly.  Primi¬ 
tive  shapes  or  volumes  (I.e.,  cubes,  boxes,  cylinders,  and  spheres)  are 
used  to  assemble  models  of  more  complicated  shapes.  These  volumes,  like 
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the  wireframes,  are  scaled  and  given  some  special  orientation  which  do 
not  necessarily  have  any  connectivity  among  them.  With  this  scheme,  the 
time  needed  for  assembly  of  a  geometric  model  is  significantly  decreased 
and  shading  for  visual  effect  is  straightforward,  since  surfaces  exist 
naturally  with  this  method. 

The  finite  element  mesh  generation  initially  doesn't  appear  to  be 
any  easier  than  with  the  wireframe  geometric  model  scheme.  The  mesh 
generation  is  restilcted  to  the  primitive  volumes.  There  is  some 
knowledge  of  the  nodes  and  elements  on  the  surface  when  working  on  the 
adjacent  volume  so  generation  of  Interface  transitions  is  easier.  The 
difficulty  lies  in  the  fact  that  transition  meshes  have  to  be  created. 

The  mesh  generation  is  constrained  at  the  common  surfaces  by  the  fact 
that  the  generation ' is  only  in  the  primitive  volumes  (l.e.,  nodes  must 
lie  on  the  surfaces).  Therefore,  the  finite  element  meshes  are  easier 
to  generate  but  still  rely  on  generating  oddly  shaped  elements  for 
transition  interfaces. 

A  third  scheme  is  referred  to  as  a  true  solids  geometric  model.  It 
has  the  knowledge  of  what  is  solid  and  what  is  vacant.  Although  the 
assembly  of  the  solid  model  may  be  by  the  use, of  primitive  volumes,  the 
common  surfaces  and,  hence,  the  common  surface  restrictions  are  eliminated., 
The  finite  element  mesh  generation  "flows"  through  the  common  surfaces 
(i.e.,  nodes  are  not  restricted  to  lie  on  the  cbmmon  surface),  thereby 
eliminating  the  use  of  transition  meshes.  The  finite  element  mesh  gemra- 
tlon  will  be  more  susceptible  to  errors  since  the  distribution  of  elements 
along  the  surfaces  have  to  be  thoTight  out  prior  to  data  input. 
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Appendix  C 


SUBSTRUCTURE  THEORY 


It  is  convenient  to  partition  a  structural  model  along  natural 
physical  subdivisions  called  substructure  models,  although  the  defini¬ 
tion  of  partitions  is  arbitrary.  For  example,  a  drydock  may  be  parti¬ 
tioned  as  shown  in  Figure  C-1.  This  figure,  illustrates  the  following 
quantities  which  are  basic  to  substructure  theory:  substructures, 
boundary  interface  nodes,  and  interior  nodes.  There  are  five  sub¬ 
structures  constituting  the-  drydock  s>.stem.  Node  b  is  a  boundary 
interface  node  betweep  contiguous  substructures  1  and  2,  and  node  i  is 
Internal  to  substructure  3  which  is  contiguous  only  with  substructure  2 
as  shown. 

Following  Przemienlecki  (1968),  we  have  the  following  vector  or 
column  matrix  definitions:' 

=  External  displacements  of  l>oundnry  interface  nodes  between 
substructures. ,  . 

Pj  =  External  displacements  of  internal  nodes  of  substructures. 

=  External  forces  on  boundary  interface  nodes. 

R.  =  External  forces  on  internal  nodes. 

The  stiffness  relations  for  the  entire  structure  expressed  in  partitioned 
matrix  form  are: 
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or,  K  D  =  R 

Substructure  analysis  Is  separated  Into  two  distinct  cases  as 
follows: 


I.  Fixed  boundary  analysis,  where  all  boundary  Interface  dof  are 
constrained  to  2ero.  The  entire  structure  is  subjected  to  external  forces. 


II.  Correction  analysis,  where'  the  boundary  Interface  dof  are 

relaxed.  The  structure  Is  subjected  to  forces  (R,  -  r/  ) ,  where  R^  is 

’  ~b  ~br  ~br 

the  vector  of  boundary  Interface  reaction  forces  at  the  artificial 
constraints,  which  are  computed  in  Case  I. 


The  total  solution  Is  the  superposition*  of  the  two  solutions: 
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*If  superposition  Is  used  within  the  context  of  Unearlxed  steps  of  a 
nonlinear  analysis  procedure,  then  substructure  theory  Is  also  valid 
for  nonlinear  analysis. 
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Considering  the  n  substructure  as  a  free  body  (i.e.,  considering  con¬ 
strained  coordinates),  and  writing  the  partitioned  substructure  stiffness 
relations,  we  get: 


(n)  '  (n) 

~bb___'__~bl_ 
(n)  I  „(n) 
~lb  I  ~il 


By  applying  the  procedure  known  as  static  condensation  to  eliminate  the 
Internal  dof  and  retaining  the  boundary  interface  dof,  we  get  the  sub¬ 
structure  stiffness  relations  (in  constrained  coordinates). 


~b  “-b  ~b 


where  the  substructure  stiffness  matrix  is: 


v(n)  _  ,f(n)  „(n)  f„(n))  „(n) 
Sb  ~  ~bb  ■  £bl  \~ii  )  ~lb 


and  reaction  forces  are: 


(s'f )  bS"’ 


The  reaction  forces  may  be  considered  a  consistent  redistribution  of 
forces  R^^^  to  the  retained  boundary  interface  nodes,  and  the  substructure 
stiffness  matrix  may  be  considered  a  super  element  stiffness  matrix. 

Now,  apply  the  Case  I  conditions  to  each  substructure,  in  turn.  That 
is,  fix  the  boundary  Interface  dof,.  thereby  isolating  substructure  n. 


r,(n)  = 

~b 


and  apply  any  external  forces  acting  on  internal  nodes 
corresponding  solution  for  internal  displacements  is: 


Also,  the  corresponding  solution  for  boundary  Interface  reactions  Is: 


5br  “  ~bi  \~il  ),~i 

These  reactions  are  necessary  to  maintain  the  constraint  condition 


The  structure  stiffness  relations  expressed  In  the  boundary  inter¬ 
face  coordinates  are  obtained  by  the  direct  stiffness  procedure  using 
the  substructure  super  elements  as  the  basic  building  blocks.  In  effect, 
a  substructure  Is  used  In  the  same  way  as  an  Individual  finite  element 
whose  Internal  dof  have  been  statically  condensed  out  (eliminated). 
Symbolically,  the  direct  stiffness  procedure  is  represented  here  for  N 
substructures  constituting  the  entire  structural  system  by  the  expres¬ 
sions:  ' 


N 


n=l 


where  is  the  resultant  external  force  on  boundary  Interface  nodes  of 
the  structure.  The  minus  sign  converts  the  substructure  reaction  forces 
to  external  forces. 

Coupling  of  substructure  stiffness  matrices  occurs  only  where 
substructures  have  common  boundaries  (contiguous).  In  the  case  of  the 
numbered  substructures  for  the  drydock  system  in  Figure  C-1,  K,  would 
have  the  form: 
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Now,  apply  the  Case  II  conditions  to  the  structlire  boundary  interface 
dof  as  follows: 


a.  Relax  the  constraints  on  the  boundary  interface  nodes,  i.e., 

D,  becomes  unknown. 

~b 

b.  Apply  the  resultant  external  force  to  the  boundary  inter¬ 
face  nodes. 

c.  Note  that  =  0  for  Case  11  since  the  interior  nodes  are 

not  accessible. 

Upon  solving  the  linear  system: 


5b  2b  = 


2b 


the  total  displacements  of  the  structure  boiindary  Interface  nodes, 

D,  ,  are  obtained.  The  values  associated  with  each  substructure, 

~b  ~bl 

can  be  picked  from  using  code  numbers  for  substructure  n  that  are 
used  in  the  direct  stiffness  assembly  procedtire. 

The  interior  node  displacements  for  each  substructure  can  be  found 
from  the  substructure  stiffness  relations: 
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Appendix  D 

EQUIVAIBNT  PIJVTE  THICKNESS  FOR  A  STEEL  CAISSON 


The  caisson  for  drydock  No.  6,  Puget  Sound  NSY,  consists  of, an 
Internal  frame  cladded  with  plates  stiffened  in  the  long  direction.  It 
was  decided  to  develop  an  equivalent  thickness  based  upon  I,  the  moment 
of  Inertia  (as  a  beam)  Instead  of  D,  the  plate  stiffness.  In  determining 
I,  the  Internal  structure  Is  Ignored  since  It  is  difficult  to  Incorporate 
the  internal  structure  In  a  rational  way.  The  stiffened  outer  shell  has 
been  designed  to  carry  the  bending  along  the  long  direction  and  is  used 
in  establishing  the  equivalent  thickness.  The  moment  of  Inertia  is  cal¬ 
culated  for  this  section. via  the  parallel  axis  theorem  using  the  centerline 
of  the  structure.  The  resulting  moments  of  inertia  are  given  In  Table 
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Table  D-1.  Moments  of  Inertia  for  Equivalent 
Plate  Thickness 


Description 

I  (In.^) 

Shell 

652,051 

Small  Angles 

90,476 

Large  Angle 

104,044 

Total 

846,571 

The  equivalent  I  Is  taken  as  the  total  and  the  thickness  Is  determined 
by  assuming  the  equivalent  solid  section  to  be  rectangular.  The  thickness 
is  calculated  as; 

h^  =  12I/b  =  196,306 

h  =  58. 1  Inches 
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